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ABSTRACT 


This  investigation  was  concerned  with  various  digital  filtering 
techniques  and  associated  constraints.  The  synthesis  procedures 
discussed  in  this  report  emphasize  the  interplay  of  the  various  critical 
design  parameters.  Generalized  design  procedures  for  lowpass,  bandpass 
and  band  stop  filters  are  developed  using  a  tabular  procedure  which 
enables  one  to  obtain  the  digital  filter  coefficients  by  inspection. 

This  approach  allows  for  a  simpler  evaluation  and  interpretation  of  such 
problems  as  coefficient  truncation,  stability  cUiu  error  constrsints  *is 
well  as  illustrating  the  importance  and  significance  of  the  concept  of 
normalization  in  digital  filters.  The  inter-relations  among  the 
foregoing  are  discussed  in  detail  leading  to  performance  curves  for 
various  implementations.  A  bandpass  and  band  stop  synthesis  technique 
which  is  accomplished  through  a  simple  conversion  of  the  lowpass 
coefficients  is  also  developed.  Bandpass  filters  having  arithmatic 
symmetry  are  then  synthesized  using  a  frequency  shift  teclinique  as  well 
as  a  lowpass  to  bandpass  transformation.  The  validation  of  these 
approaches  for  various  ratios  of  sampling  rate  to  carrier  frequency  is 
discussed.  An  analysis  of  synthesis  errors  is  then  accomplished.  Under 
the  assumption  that  tabular  data  is  available,  design  procedures  which 
minimize  the  sum- squared  error  are  developed  for  design  of  non- recursive 
digital  filters.  A  second  approach  to  the  design  of  these  filters 
was  accomplished  under  the  assumption  that  a  satisfactory  recursive 
digital  filter  design  using  the  bilinear  transform  was  available. 


iii 


TABLE  OF  CONTENTS 


Page 

INTRODUCTION  1 

A.  Summary  of  Results  2 

GENERAL  DESIGN  PROCEDURES  7 

A.  Bilinear  Transformation  -  Lowpass  Applications  7 

B.  Bandpass  Digital  Filters  -  LP-BP  Transformation  18 

C.  Bandpass  Digital  Filters  -  A  Shifting  Technique  25 

D.  Bandstop  Digital  Filters  -  A  Feedback  Technique  28 

ARITHMETICALLY  SYMMETRICAL  BANDPASS  FILTERS  34 

A.  Design  Procedure  for  Elliptic  Analog  Filters  34 

B.  Digitized  Lowpass  Elliptic  Filters  37 

C.  Application  of  the  Shifting  Technique  44 

D.  Symmetry  Errors  and  Ch'  me  of  Center  Frequency  47 

E.  Application  of  the  LP-BP  Transformation  53 

NONRECURSIVE  DIGITAL  FILTERS  60 

A.  Nonrecursive  Design  From  Tabulated  Data  62 

B.  A  Digital  Impulse  Invariant  Technique  65 

C.  Zero  Removal  and  Relocation  7  1 

ERROR  CONSIDERATIONS  78 

A.  Stability  and  Coefficient  Accuracy  78 

B.  Cascade  vs.  Direct  Synthesis  84 

C.  Computational  Quantization  93 

D.  An  Auxiliary  Storage  Technique  98 

RECOMMENDATIONS  108 


v 


110 


APPENDIX  A  A  BANDPASS  SAMPLING  TECHNIQUE 
APPENDIX  B  DIGITAL  OSCILLATORS 
APPENDIX  C  REFERENCES 
APPENDIX  D  BIBLIOGRAPHY 


EVALUATION 


The  increasing  speed,  as  well  as  the  decreasing  size  and  cost 
associated  with  digital  circuitry,  that  have  been  and  still  are 
resulting  from  the  advances  in  the  field  of  microelectronics,  brings 
within  the  likelihood  that  digital  filters  will  perform  within 

real  time  devices,  almost  all  the  functions  now  performed  with  analog 
components.  The  advantage  of  these  digital  filters  in  terms  of  the 
increased  accuracies  attainable,  the  ability  to  change  filter  shapes, 
the  time  sharing  capabilities,  the  simplicity  of  the  components  needed 
(multipliers,  adders  and  storage  devices)  and  the  variety  of  attainable 
filter  shapes,  among  many  other  advantages,  all  serve  to  illustrate  the 
importance  and  widespread  application  of  this  discipline. 

Apart  from  "The  Advanced  Digital  Processing"  Contract  now  under  way 
(Contract  No.  F30602-69-C-0199)  it  is  strongly  recoimended  that  further 
work  exploit  the  rich  opportunities  of  this  fruitful  area.  In  particular 
the  error  reduction  technique  discussed  in  Section  V-D  has  shown 
sufficient  merit  to  warrant  experimental  breadboarding.  Furthermore 
an  experimental  investigation  should  be  initiated  concerning  the 
digital  oscillator  implementation  technique  discussed  in  Appendix  B. 

A  further  area  in  this  program  looks  highly  rewarding  at  this  time, 
and  that  is  the  question  of  direct  digital  design  in  complement  to 
the  presently  developing  family  of  digital  equivalence  procedures. 

This  approach  (direct  digital  synthesis)  warrants  detailed 
investigation  in  that  it  will  allow  for  more  flexibility  in  the  designs 
that  can  be  achieved. 

CHARLES  N.  MEYER 
Project  Engineer 
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I.  INTRODUCTION 


The  field  of  digital  filters  is  based  essentially  on  the  mathematics  of 
difference  equations.  In  addition,  the  work  accomplished  in  the  fields  of  sam¬ 
pling,  quantizing  and  related  transformation  techniques,  have  helped  in  both 
the  reinterpretation  and  furthering  of  this  discipline. 

The  increasing  speed,  as  well  as  the  decreasing  size  and  cost  associat¬ 
ed  with  digital  circuitry,  that  have  been  and  still  are  resulting  from  the 
advances  in  the  field  of  microelectronics,  brings  within  view  the  likelihood  that 
digita'  filters  will  oerform  within  real-time  devices,  almost  all  the  functions 
now  performed  with  analog  components.  The  advantages  of  these  digital  filters 
in  terms  of  the  increased  accuracies  attainable,  the  ability  to  change  filter 
shapes,  the  time  sharing  capabilities,  the  simplicity  of  the  components  needed 
(multipliers,  adders  and  storage  devices)  and  the  variety  of  attainable  filter 
shapco,  among  many  other  advantages,  all  serve  to  illustrate  the  importance 
and  widespread  application  of  this  discipline. 

Thus,  it  is  important  that  design  procedures  for  digital  filters  be  devel¬ 
oped  so  that  filters  can  be  practically  and  economically  syn'hesized,  and  that 
eng ineer ing  designs  can  be  accomplished. 

This  investigation  v/as  concerned  with  various  digital  filtering  techniques 
and  associated  constraints.  Section  II  through  Section  V  develop  in  detail  the 
foregoing.  Recommendations  for  further  study  and  development  are  discussed 
in  Section  VI.  The  essential  results  of  this  investigation  are  summarized 
be  1  ow . 
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A.  SUMMARY  OF  RESULTS 


Section  II  -  General  Design  Procedures 

Design  procedures  for  lowpass,  bandpass  and  bandstop  digital  filters, 
are  discussed.  The  bilinear  transformation  is  applied  to  the  design  of  lowpass 
fil'.ers  in  such  a  way  that  the  relationships  among  the  required  sampling  rate, 
upper  cutoff  frequency  and  order  of  the  filter  are  brought  out.  The  approach 
illustrates  the  importance  of  the  ratio  of  sampling  to  cutoff  frequency.  The  con¬ 
cept  of  normalization  and  its  significance  is  also  discussed.  A  tabular  proce¬ 
dure  is  developed  where  these  digital  coefficients  can  be  written  down  by 
inspection.  Furthermore,  this  approach  allows  for  a  simpler  evaluation  and 
interpretation  of  coefficient  truncation,  stability  and  errors,  as  discussed  in 
detail  in  Section  V. 

A  lowpass  to  bandpass  transformation  is  used  in  conjunction  with  the 
bilinear  transformation.  This  combined  transformation  is  then  structured  in 
such  a  form  that  a  bandpass  design  can  be  accomplished  by  spec  lfying  the  ratio 
of  sampling  frequency  to  bandwidth  and  sampling  frequency  to  carrier  frequency. 
Thus,  once  again,  normalized  frequency  parameters  are  shown  to  be  significant. 
An  example  illuo  crates  how  the  foregoing  are  related  to  the  factor  Q,  an  often 
used  parameter  in  analog  filters.  By  an  adjustment  of  the  sampling -to-carrier 
frequency  ratio  a  technique  is  described  which  allows  a  lowpass  filter  to  be 
converted  to  a  bandpass  design  without  changing  any  of  the  coefficients. 

A  second  bandpass  design  procedure  which  utilizes  a  frequency  shifting 
technique  is  described.  This  approach  shows  considerable  promise  in  the  design 
of  arithmetically  symmetrical  filters  as  discussed  in  Section  III.  This  shifting 
procedure  also  suggests  a  technique  for  generating  the  coefficients  "on  line" 
from  a  lowpass  design. 

Lastly,  a  feedback  technique  is  described  which  offers  a  simple 
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unce  again,  the  procedure 


procedure  for  designing  bandstop  or  notch  fillers, 
is  accomplished  through  a  simple  conversion  of  the  lowpa68  coefficients  and 
shows  clearly  the  relationship  between  the  depth  of  the  notch  and  value  of  the 
resulting  coefficients. 

Section  III  -  Arithmetically  Symmetrical  Bandpass  Filters 

A  formalized,  step-by-step  procedure  is  outlined  for  the  design  of 
elliptic  analog  filters  with  examples  illustrating  the  technique.  These  filters 
are  then  digitized  as  lowpass  filters  using  the  approach  discussed  in  Section  II. 
To  convert  these  digital  lowpass  filters  to  bandpass,  the  frequency  shifting 
technique  of  Section  II  is  applied  and  the  mechanism  of  the  coefficient  changes 
described.  The  examples  chosen  indicated  that  it  was  possible  to  reduce  the 
symmetry  errors  to  zero.  Thus,  a  mathematical  analysis  of  symmetry  errors 
utilizing  this  shifting  technique  was  accomplished.  It  was  shown  that  there  are 
two  types  of  errors.  A  symmetry  error  about  the  carrier  frequency  and  the 
error  resulting  from  the  shift  from  lowpass  to  bandpass.  It  is  shown  that  the 
former  error  can  be  reduced  to  zero  and  the  latter  error  to  a  negligible  amount 
by  appropriately  choosing  the  relationship  between  carrier  frequency  and  sam¬ 
pling  rate. 

It  is  then  shown  that  for  a  particular  ratio  of  sampling  to  carrier  fre¬ 
quency,  the  I.P-BP  transformation  yielded  an  arithmetically  symmetrical 
design  with  zero  symmetry  error.  This  technique  is  then  compared  with  the 
shifting  technique  at  the  same  ratio.  Although  in  general  the  shifting  technique 
is  excellently  suited  for  the  design  of  symmetrical  filters,  (and  the  LP-BP 
te  chnique  is  not)  the  LP-BP  transformation  at  tliis  one  sampling -to -car rier 
frequency  ratio  proves  to  be  a  simpler  technique  to  apply  if  one  alread'  has  the 
digital  lowpass  filter  implementation. 
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Section  IV  -  Nonrecursive  Digital  Filters 

Under  the  assumption  that  tabular  data  is  available  in  terms  of  the 
desired  frequency  characteristics,  (both  amplitude  and  phase),  a  design  pro¬ 
cedure  for  nonrecursive  digital  filters  is  developed  which  minimizes  the  sum 
squared  error  between  the  desired  and  actual  response  characteristics.  These 
design  equations  are  then  related  to  an  equivalent  procedure  which  can  be  used 
when  the  desired  response  characteristics  are  available  as  a  continuous  function 
of  frequency. 

A  second  approach  to  the  design  of  nonrecursive  filters  was  accomplished 
under  the  assumption  that  a  satisfactory  recursive  digital  filter  design  was 
available. 

Thus  the  procedure  consists  essentially  of  developing  a  finite  Fourier 
series  approximation  to  the  recursive  filter.  It  was  assumed  that  the  recursive 
filter  was  designed  using  the  bilinear  transform  and  that  the  nonrecursive 
approximation  was  to  be  a  minimum  mean  square  approximation  to  the  recursive 
form.  A  simple  technique  (referred  to  as  a  digital  impulse  invariant  technique) 
for  converting  the  recursive  coefficients  to  the  required  nonrecursive  coefficients 
is  developed.  This  approach  yields  zero  error  in  the  impulse  responses  of  the 
two  filters  up  tc  the  number  of  terms  retained  in  the  nonrecursive  form. 

It  was  noted  that  there  are  essentially  two  types  of  errors  in  nonrecur¬ 
sive  filter  design.  The  first  was  due  to  a  transformation  procedure,  the  second 
due  to  truncation.  The  process  of  zero  removal  and  relocation  is  discussed  in 
the  context  of  removing  this  second  source  of  error.  The  interpretation  of 
raised  cosine  pulses  is  accomplished  using  this  technique  a£  well  as  the  relation¬ 
ships  among  several  other  nonrecursive  design  procedures.  Ideally,  this 
procedure  offers  the  possibility  of  a  complete  characterization  of  finite  pulses. 
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Section  V  -  Error  Considerations 


It  is  noted  that  there  are  three  types  of  errors  which  occur  in  digital 
filter  processing.  The  first  of  these  is  related  to  the  process  of  sampling  and 
quantizing  the  input  signal.  The  second  occurs  as  a  result  of  truncating  the 
representation  of  the  filter  coefficients.  The  third  error  is  referred  to  as 
computational  quantization.  That  is,  the  errors  that  occur  as  a  result  of 
quantizing  the  weighced  multiplications  and  additions  that  occur  within  the  digital 
filter's  arithmetic  unit.  In  recursive  filters  the  results  of  these  computations 
are  fed  back  and  utilized  in  later  computations.  This  section  discusses  the 
latter  two  types  of  errors. 

The  relationship  between  coefficient  accuracy  and  filter  stability  is 
determined  using  the  tabular  approach  discussed  in  Section  II.  Simple  relation¬ 
ships  are  derived  which  clearly  show  the  interplay  between  the  order  of  the 
filter,  the  ratio  of  sampling  rate  to  critical  filter  parameters  and  the  bit 
requirements  for  the  filter's  coefficients.  Using  this  tabular  approach,  the 
filter  degradation  and  onset  of  instability  is  interpreted  as  a  function  of  the  bit 
truncation  of  the  coefficients.  It  is  note  .nat  when  the  ratio  of  sampling  rate 
to  filter  cutoff  frequency  (in  conjunction  with  the  order  of  the  filter)  go  beyond 
a  specified  value  relative  to  the  number  of  bits  retained,  the  direct  approach 
should  not  be  used. 

The  foregoing  constraints  are  illustrated  through  the  design  of  a  high 
order  lowpass  filter  synthesized  using  the  bilinear  transformation.  The  supeii- 
ority  of  the  cascade  approach  (and  the  conditions  under  which  it  is  superior)  is 
then  illustrated  by  synthesizing  a  filter  i  tilizing  the  same  number  of  bits  per 
coefficient  as  in  the  direct  form.  This  is  accomplished  twice.  The  first  time 
the  critical  parameter  are  adjusted  so  that  the  direct  and  cascade  approach 
yield  similar  results.  Then  this  same  is  parameter  adjusted  so  that  the  direct 
form  docs  not  yield  useful  results  whereas  the  cascade  approach  does. 
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The  relationship  between  the  direct  and  canonical  recursion  implemen¬ 
tations  are  discussed  as  to  their  performance  when  computational  quantization 
occurs . 

Illustrative  examples  indicate  that  a  crossover  point  exists  at  which  one 
implementation  is  more  preferable  than  the  other.  It  is  noted  once  again  that 
the  determining  factor  is  the  normalized  sampling  rate. 

Lastly,  an  auxiliary  storage  technique  is  discussed  which  offers  the 
possibility  of  reducing  the  errors  due  to  computational  quantization.  An  exam¬ 
ple  illustrates  the  approach. 

Appendix  A  -  A  Bandpass  Sampling  Technique 

Allowable  sampling  rates  are  discussed  for  bandpass  signals  so  that  no 
overlap  distortion  occurs.  A  curve  is  obtained  relating  the  allowable  sampling 
rate  to  the  bandwidth  and  carrier  frequency.  As  noted  in  the  previous  sections, 
these  relationships  are  of  extreme  importance  for  bandpass  filter  design. 

Appendix  B  -  Digital  Oscillators 

Various  implementations  for  a  digital  oscillator  are  discussed.  An 
implementation  utilizing  the  canonical  form  representation  yields  the  smallest 
errors  (due  to  computational  quantization)  and  also  requires  the  least  amount 
of  hardware.  This  device  is  exceedingly  useful  both  as  a  frequency  synthesizer 
and  as  a  device  to  be  used  in  conjunction  with  tracking  filters,  frequency 
translators  and  bandpass  filtering. 
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II.  GENERAL  DESIGN  PROCEDURES 


The  general  problem  of  recursive  digital  filter  design  is  the  determina¬ 
tion  of  the  coefficients  a.  and  b.  so  that  a  desired  filter  characteristic  is  obtained. 

J  J 

A  basic  approach  to  this  problem  has  been  to  utilize  existing  continuous  (analog) 

filter  theory  to  "find"  a  suitable  response  and  then  apply  a  transformation 

technique  which  digitizes  the  continuous  filter  in  such  a  waj  that  the  response 

& 

is  transformed  without  any  appreciable  distortion.  This  is  referred  to  as  a 
digital  equivalence  technique. 

A  transformation  which  has  been  used  with  considerable  success  is  the 
so-called  bilinear  ^  or  Tustin^^  transformation  given  by  the  mapping  between 
the  variables  S  and  Z  *  given  by 


S  — ■ 


_2 

T 


(1  -  z~j) 

(i  +  z"1) 


II-  3 


or 


„-l_  1  -  (ST/2) 

z  i  +  (ST/zy 


II  -  4 


The  advantageous  properties  of  this  transformation  ir.  terms  of  preser¬ 
vation  of  stability  (the  left  of  the  s  plane  is  mapped  into  the  interior  of  the  unit 
circle  of  the  Z  plane),  the  maintenance  of  the  cascading  property  (as  well  as 
dc  gain)  and  the  simplicity  of  its  application  in  that  it  is  purely  algebraic  in 
nature,  has  been  discussed  by  several  authors.' ^  A  disadvantage  of  this 
technique  is  that  the  transformation  given  in  II - 3  causes  a  distortion  in  the 
frequency  domain.  That  is,  with  s  =  j  <o  and  Z  *  =  e  II _ 3  becomes 


wT 


— »  tan  to  T/2 


II -5 


The  approach  for  the  design  of  nonrecursive  filters  often  uses  a  more 
direct  approach  without  resorting  to  transformation  techniques.  These  filters 
(b.  =  o,i<  j  <  M  )  will  be  discussed  in  Section  IV. 
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Thus,  the  relationship  between  the  analog  and  digital  frequency  variables  is  not 
linear.  The  deviation  from  linearity  given  by  II- 5  is  dependent  on  the  pr oduct 
(u>T/2).  Since  T  is  the  reciprocal  of  the  sampling  rate,  this  distortion  or  so- 
called  warping  of  the  frequency  scale  is  reduced  when  the  ratio  of  sampling  rate 
to  critical  frequency  points  is  made  larger.  There  are  however,  other  con¬ 
straints  on  this  ratio  which  will  be  discussed.  Furthermore,  compensation  of 
this  distortion  will  also  be  discussed. 

In  that  the  bilinear  transformation  properly  occupies  a  central  position 
in  the  design  of  digital  filters,  it  has  been  utilized  extensively  in  this  investiga¬ 
tion.  In  order  to  illustrate  the  design  approach  and  discuss  the  interrelation¬ 
ships  among  the  various  critical  parameters,  consider  a  normalized  (unity 
cutoff  frequency,  in  rad/sec)  low  pass  (LP)  analog  (or  continuous)  filter  of  the 
form 

N 

V  A  Sn 

lj  n 

C(S)  =  Ji^° -  11-6* 

y  B  s" 

-j  n 
n  =  0 

The  concept  of  a  normalized  frequency  response  characteristic  is  widely  accept¬ 
ed  in  the  synthesis  of  continuous  filters  as  a  procedure  which  allows  for  a 
universal  design.  In  order  to  convert  this  normalized  response  to  a  low-pass 
filter  with  an  upper  cutoff  frequency  of  oj  ,  if  one  applies  the  conventional  LP  to 
LP  transformation 


Although  the  numerator  and  denominator  are  usually  ol  -'iflerent  order, 
they  can  always  be  written  as  shown,  by  adding  the  required  zeros. 


in  conjunction  with  the  bilinear  transformation  from  S  to  Z,  it  can  be  shown 
that  the  synthesized  digital  filter  H(Z)  becomes 


H(Z)  =  G(A) 

U) 

u 


i  -  z 


1  +  z 


n-8 


or 


H(Z)  = 


N  T 

£  An(^-)N-n  (1  -  Z'1)11  (1  +  Z^)^" 
n  =  0 

N  w  T 

£  Bn(-y-)N-n  (1  -  Z_1)n  (1  +  Z_1)N'n 
n  =  0 


II -9 


In  the  determination  of  II-9,  fractions  were  cleared  by  multiplying  the  numerator 

.  N  -IN  — 

and  denominator  by  (ui^T/2)  and  (1  +  Z  )  .  Now,  suppose  is  the  desired 
digital  upper  cutoff  frequency  of  the  digital  filter. 

Using  the  transformation  of  II-5,  the  analog  cutoff  frequency  variable 
is  replaced  by  its  equivalence  in  terms  of  and  then  II - 9  becomes 


N  N 

N-r.  ,,  ,-l,n  „  .  _-l,N-n 


H(Z)  - 


V  a  Z~n 

— /  n 
n  =  0 


V  A  K^'H  (1  -  Z_,)n  (1  +  Z'1) 
i-J  n  ui 

n  =  0  u 

_T_ - -  TJ - 

y  B  K^‘n  (1  -  Z'1)”  (1  +  Z"])N'n  v  b  z'n 
Li  n  i-i  n 

n= 0  U  n  = 0 


II  - 1  0 


where  the  a  and  b  are  the  coefficients  of  Z  n  obtain;1'!  bv  expanding  the  sum- 
n  n 

mationa  in  II- 10  and 


—  T 

K*  —  tan(  w  ■  w  ) 
u>  u  2 


II- 11 


Thus,  the  frequency  response  of  a  digital  filter  is  a  function  of  the  normalized 
(with  respect  to  the  sampling  rate)  digital  cutoff  frequency.  Furthermore,  since 
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Z  *  =  e  and  the  A  and  B  are  the  constants  obtained  from  the  normalized 

n  n 

lowpaas  analog  filter,  the  entire  digital  frequency  characteristic  (H(Z))  is  such 
that  the  digital  frequency  variable  uappears  only  as  the  product  uT.  Thus  if  this 
characteristic  is  plotted  as  a  function  of  u  T  and  not  u,  one  obtains  a  normalized 
digital  frequency  characteristic  which  provides  a  universal  curve  in  the  same 
sense  as  normalization  with  respect  to  u  =  1  rad/sec  did  in  the  analog  case. 

This  is  exceedingly  important  to  note  in  that 


the  coefficients  (an,  bn)  for  this  digital  lowpass 

filter  remain  fixed  as  long  as  the  product  w  T  = 

=  2  it  r  /f  remains  fi'-cd.  u 

u  s 

Therefore,  the  coefficients  for  a  digital  lowpass  filter  with  a  cutoff  at 
7^  =  10  KHz  operating  on  data  sampled  at  fg  =  100  KHz  are  exactly  the  same  as 
those  obtained  for  any  other  T  and  f  as  long  as  the  ratio  T/f  =  ,  1.  Hence, 
by  plotting  the  characteristic  as  a  function  of  w  T  (i.e.  in  radians  or  degrees) 
one  has  a  "universal"  digital  filter  characteristic  (the  digital  coefficients  re¬ 
main  fixed)  for  that  particular  ratio  of T  to  f  (or  equivalently  a  fixed  K—  ). 

'‘’u 

To  convert  this  angular  abscissa  to  actual  frequencies,  one  merely  utilizes  the 
actual  sampling  rates  used  on  the  input  data.  These  rates  are  of  course  dictated 
by  the  bandwidth  of  the  input  data  and  the  constraints  of  the  sampling  theorem. 
For  the  foregoing  reasons  all  digital  filter  characteristics  in  this  report  arc 
plotted  as  a  function  of  angle.  The  above  results  are  of  course  applicable  to 
bandpass,  bandstop  and  all  other  digital  filter  forms.  The  critical  ratio  and 
parameters  in  these  cases  v,ill  be  discussed  in  II-B. 

In  order  to  obtain  the  a's  and  b's  of  the  digital  filter  in  terms  of  K—  and 

u 

the  A's  and  B's  of  the  analog  filter,  11-10  must  be  expanded.  It  is  also  to  be 
noted  that  the  form  of  the  numerator  and  denominator  of  II- 10  is  similar  to  one 
another.  Thus  the  expansion  equations  will  likewise  be  similar.  During  this 
investigation  it  has  been  shown  that  these  digital  coefficients  for  any  order  filter 
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f 


are  given  by  the  following  table. 


TABLE  I 


Coefficients 

P  K.- 

U  u) 

u 

p,^-1 

u 

P.K1!-2 

2  w 

u 

.  .  .  . 

PN 

K  — 

L) 

U 

z0<aobo} 

i 

i 

i 

.  .  .  . 

l 

(  ) 

(  ) 

( 

) 

-2 

z 

<=? 

<  > 

<  ) 

( 

) 

•  • 

• 

• 

d. 

1  •  J 

Z'N<aN*V 

(  ) 

( 

) 

where 


CN  = 


N 


r  r  (N  -  r)  1 
The  entry  cL  ^  in  row  i  and  column  j  is  given  by 

d.  .  =  d.  -  fd  +  d  .  ] 

i .  .1  i  -  j  - 1  L  1 "  1  -  J  -  1  i  -  1 .  J  J 


!I  -  1 2 


The  interpretation  of  this  table  to  obtain  the  digital  filter  coefficients  is  as 
follows:  To  obtain  the  numerator  coefficients  a,  ,  substitute  for  P.  the  appro- 

K.  J 

priately  subscripted  numerator  coefficient  of  the  analog  filter  and  sunt  all  the 
products  of  the  table  entries  and  column  headings.  The  same  .procedure  is  used 
for  the  denominator  coefficients  except  that  the  analog  denominator  coefficients 
are  now  substituted  for  the  F  ^ . 

As  noted  from  this  Tabic,  the  first  row  consists  of  all  "ones11  and  the 
first  column  consists  of  the  binomial  coefficients.  It  can  also  be  shown  that 
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the  last  row  ie  identical  to  the  first  row  except  with  alternating  signs  and  the 
last  column  identical  to  the  first  column  except  for  alternating  signs.  Utilizing 
11-12,  it  can  also  be  shown  that  the  sum  of  the  table  entries  over  any  column 
except  the  first  is  equal  to  Zero.  The  above  relationships  \-<ill  be  found  to  be 
useful  in  the  discussion  of  coefficient  accuracy  and  stability  considerations  dis¬ 
cussed  in  Section  V. 

As  an  example  of  the  above,  Table  II  lists  the  digital  filter  coefficients 
for  a  sixth  order  digital  filter.  The  coefficients  in  this  table  are  obtained 
utilizing  II  -  1 2  . 


TABLE  II 


Sixth  Order  Digital  Filter  Coefficients 


-  2  P.KL  -  4  P„K—  -  6P, 
4  ui  o  ui  6 

u  u 

(a,  or  b,)  =  15  P_K^L  +  5  P.k!  -  P,k1  -  3  P.kL  -  P.kL+5P,K-  +  1  5  P, 
Z  Z  0  uj  1  w  Z  w  3w  4  to  5  w  t> 

u  u  u  u  u  u 


(a  or  bj  =20P.K- 

5  5  U  u) 


-  4  P,KL 


+  4P.K- 

4  ui 


-  20  P 


(a,  or  b.)  =  13P„K6  -  5P.K--  P,k1  3  P,k1  -  P.K1-5P..K-  +  15P 

4  4  Ocj  1  u>  3oj  4  u>  5  w 

U  \1  U  U  U 


(ac  or  b,)  =  6  PnK-  -  -1  P.k!+  2  P,Ki 

D  0  U  u)  iU>  <-  W 

u  u  u 


-2P/1+4P.K-  -  6  P/ 
4  u;  5  W  6 

u  u 


(a ,  or  b, )  =  PnK—  -  P,k!+  P9k1  *  P,kL  +  P.kI-  PtK-  +  P, 

0  6  0  w  1  2uj  3  ^  4  u  5  w  6 

u  u  u  u  u  u 


Although  error  considerations  will  be  discussed  in  detail  in  Section  V, 
Table  II  illustrates  some  ot  the  constraints  among  the  various  parameters  and 
how  they  affect  the  required  number  of  bits  to  be  retained  for  the  coefficients. 
As  can  be  seen  from  this  Table  (and  Table  I),  K—  appears  in  the  equations 
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raised  to  powers  up  to  the  order  of  the  filter.  This  factor  (given  by  II - 1  1 )  is 

proportional  to  the  ratio  of  upper  cutoff  frequency  to  sampling  frequency.  Thus, 

N 

when  this  ratio  is  small  and  the  order  of  the  filter  high,  K—  is  a  small  number. 

u 

Qualitatively  from  these  tables  it  can  be  seen  that  to  include  the  effect  of  Pq  a 

sufficient  number  of  bits  must  be  retained.  For  example,  if  the  ratio  of 

sampling  to  cutoff  frequency  is  30:1  then  K—  ax  1  0  ^  and  as  many  as  six  decimal 

wu 

digits  might  be  required  in  the  synthesis,  It  should  be  recalled  mat  the  sampling 
rate  is  dictated  by  the  signal  bandwidth  and  aliasing  considerations  and  there¬ 
fore  cannot  be  made  arbitrarily  smaller.  Thus,  if  it  is  required  to  design  a 
filter  with  a  3  db  point  at  30  KH*  to  filter  a  signal  with  a  bandwidth  of  .  3  MHt, 
then  with  a  sampling  rate  of  3  times  maximum  the  ratio  discussed  above  becomes 
30.1  =  f  :T  . 

8  U 

As  an  illustrative  example  of  the  use  of  the  foregoing  tables,  consider 
the  design  of  a  normalized  lowpass  elliptic  filter  designed  to  meet  the  follow¬ 
ing  specifications. 


Passband  Spec  .  ,  01  db  ripple 


0  <  w  <  1 


Stopband  Spec.  40  db  attenuation  1  38  <  u  <  oo 


The  above  specifications  can  be  shown  to  yield  a  sixth  order  analog  filter  of  the 
form 


S4  +  A  S1^  +  A 

G  (S)  =  —5 - 5 - 3 - 3 - 2 - 

S  +  B,S  +  B,S  +  B,S5  +  B,S  +  B.S  +  Bn 
5  4  3  Z  10 


11-13 


These  filters  will  be  discussed  at  length  in  Section  III. 

Note  the  normalizations  on  the  analog  filter.  That  is  the  stopband 
begins  38%  beyond  the  passband, 
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where  the  coefficients  are  given  as 

TABLE  III 

AQ  =  6.  7  9609  B2  =  5.  7  1737 

A  .  =  5.  401  OK  B  .  =  6.  .39689 

2  :> 

B,  =  1 . 26743  B.  =  5.  06663 

0  4 

Bj  =  3. 44669  B5  =  2. 62)93 

B6=  1 

Solving  for  the  roots  of  the  numerator  and  denominator  yield  the  pole -zero  con¬ 
figuration  for  this  filter  given  by 

TABLE  IV 
z.  =  t  j  1.84538 

(zeros) 

z  =  t  j  1.41268 

p  =  10102  t  i  1.  12369 

p2  =  -.  39699  t  .9921  1  (poles) 

p3  =  -.81294  t  • 43946 

It  was  decided  to  digitize  this  filter  for  a  ratio  of  sampling  rate  to  cut¬ 
off  frequency  of  18:1.  This  could  correspond  to  the  problem  of  filtering  sampled 
data  in  the  band  from  0  to  1  KHz  from  a  signal  with  a  bandwidth  of  6  KHz  that 
has  been  sampled  at  three  times  its  bandwidth  or  f  =  18  KHz,  thus  yielding 

f  U  t.  18. 
s'  u 


The  procedure  to  obtain  these  coefficients  will  be  discussed  in  Section 
III,  where  additional  references  on  analog  filter  coefficient  determination  will 
be  given. 
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For  this  ratio,  the  product  u^T  becomes 

r 

c^T  2  IT  (rad)  =  20°  II- 14 

s 

Thus  K—  becomes 

OJ 

u 

d>  T 

K-  =■  tan(-y-  )  =  tan{  10°)  =.  17633  II- )  5 

u 

IJ9ing  the  analog  coefficients  given  in  Table  III,  the  factor  K—  given 

u 

ibove,  and  the  transformation  equations  of  Table  II,  the  digital  filter  coef¬ 
ficients  obtained  by  summing  each  row  of  thi9  table  are  given  as 

TABLE  V 

aQ  =  .0365179  b0 =  1.6605333 

aj  =  -.0505161  bj  =  -  8. 1507334 

a2=  -.0332495  b2  =  17.0485046 

a3  =  .1075630  =  - 1  9 . 39 1 2 1 46 

a  =-.0332495  b.  =  12.6265180 

4  4 

a5  =  -.0505161  b5  =  '  4- 4568329 

a,  -  .0365179  b,  =  .6656633 

o  o 

(numerator)  (denominator) 

A  plot  of  the  resulting  magnitude  characteris tie  is  shown  Figure  II - 1  . 
Two  significant  items  are  to  be  noted  from  the  above.  First,  the  digital 
filter  obtained  is  valid  (without  changing  any  coefficient)  at  any  frequency  as 
long  as  the  18  ;1  ratio  is  maintained.  Second,  the  plotting  of  the  magnitude 
characteristic  is  shown  as  a  fvnetion  of  angle  --  not  frequency.  Thus,  the  cut¬ 
off's  are  adjusted  to  the  desired  frequency  by  a  selection  of  the  sampling  rate. 

1  6 


MA.6 
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The  foregoing  illustrated  a  design  procedure  for  lowpass  filters  which 
is  simple  to  apply  and  which  yields  accurate  results.  Consideration  on  coef¬ 
ficient  bit  lengths,  errors  and  implementation  procedures,  will  be  discussed 
in  Section  V. 

B.  BANDPASS  DIGITAL  FILTERS  -  LP-BP  TRANSFORMATION 

As  noted  in  the  previous  section,  digital  filter  normalization  can  be 
accomplished  for  all  filter  forms.  A  transformation  from  low  pass  to  band¬ 
pass  will  now  be  discussed  which  illustrates  the  manner  in  which  the  critical 
parameters  affect  the  design. 

The  conventional  lowpass  to  bandpass  transformation  (LP-BP)  used  in 
analog  filter  design  is  given  as 

LP-BP  (analog) 


S  — 


■>  2 
S“  +  w 

o 


TT5" 


11-16 


where  B  =  (bandwidth),  (center  frequency  squared)  and  u>  ^  and 

are  respectively  the  lower  and  upper  cutoff  analog  frequencies.  If  11-16  is 
combined  with  the  bilinear  transformation  given  in  II - 3 ,  then  II- 16  becomes 


S  — 


1+u4T2/4 

O 

— ET - 


w2TZ  w2T2  \  .  2 

— T—  ) /O  +  -T“)  z +  z 


1 


11-17 


Now,  the  bilinear  transformation  transforms  analog  frequencies  to  digital  fre¬ 
quencies  as 


tan 


uT 


11-18 


Thus,  utilizing  the  above  and  the  relationships  between  u>q,  and  the  factors 
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in  II-  17  become 


1  +io2T2/4 
o 

— ST - 

~r 


ui  T  u>?T 

1  +  tan( — > — )  tan( — « — ) 


1 


T  u.T 

tan( — « — )  -  tan( — « — )  tan 


w.  -w  )  T 

1  ") 


where 


2  2  , 

1  -u  T  /  4 
o 

1  +  oj2T2/4 
o  ' 


1  -  tan2 


1  +  tan2 


=  cos  u  T 
o 


Thus  II -  1 7  becomes 

LP-BP  (digital) 

1  -  2(cos  u  T)Z'*  +  Z'2 

S - 2 - 72 - 

K_<1  -  Z 
B 

The  transformation  from  S  to  Z  is  now  given  in  terms  of  the  digital 
of  (1)  digital  carrier  frequency  and  (2)  digital  bandwidth  (u^  - 
again,  the  parameters  are  normalized  with  respect  to  the  sampling 


--■h  n->’ 

B 


11-20 


11-21 

parameters 
,  Once 
rate . 
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Similarly,  the  lowpass  to  bandstop  transformation  is  merely  the  inverse 


of  11-21  or 


Kb(1  -  Z'2) 

)  -  2<cos  u  T)Z'1  +  Z*2 
o 


11-22 


Thus,  one  now  has  two  simple  transformations  to  convert  a  normalized  lowpass 
analog  filter  to  either  a  bandpass  or  bandstop  digital  filter  by  specifying  the 
normalized  digital  center  frequency  and  normalized  digital  bandwidth.  Qualita¬ 
tively,  this  equation  can  be  related  to  the  lowpass  transformation  by  noting  that 

the  factor  K=  has  replaced  the  loy/pass  factor  K—  and  the  numerator  of  11-21 
B  “u 

can  be  interpreted  by  noting  that 


l  ?  l  -  j  <a  T 

1  -  2(cos  w  T)Z_i  +  Z~C  =  (Z"1  -  e  °  )<Z 


•1  +juoT 


11-23 


- 1  “j  %T  _  +  — 

Thus  since  Z  =  e  ,  the  roots  have  been  shifted  to  w  -  -  the  digital 

carrier  frequency. 

As  a  simple  example  of  the  application  of  the  LP-BS  transformation, 
consider  the  lowpass  R-C  filter 

1 


LP 


(S)  = 


1T7T 


11-24 


substituting  11-22  in  the  above  yields 

1  -  2(cos  w  T)  Z'1  +  Z'2 

H(Z)  =  - 2 - j - T  11-25 

(1  +  K  —  )  -  2(cos  WqT)  Z’1  +  (1  -  K-g-)  Z 

Thus  11-25  represents  a  digitized  R-C  bandstop  filter  for  any  normalized  digital 
center  frequency  and  normalized  digital  bandwidth  (or  notch  or  rejection  band). 
The  resulting  digital  filter  is  of  the  same  form  as  would  result  from  applying 
the  bilinear  transform  directly  to  the  familiar  analog  notch  filter  given  by 
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11-26 


~l - 1 

S  +  2  4  XS  +  X 

where  X  is  the  analog  notch  frequency  and  |  is  the  analog  damping  ratio.  This 
is  so  because  the  R-C  filter  is  converted  to  the  above  by  the  LP-BS  (analog) 
transformation. 

A  further  interpretation  of  the  factor  Kg-  or  terms  of  the  Q  of  the  filter 
is  as  follows 


Kg  =  tan  ) T / 2 J  -  tanir  B^T 

Ff  =  dig  1  bW 

Now  (l/T)  =  fg  =  sampling  rate.  However,  if  one  assumes  that 

f  =  rT 
s  o 

and 


„  o  carrier  frequency 
U  =  ^  Bandwidth 

Bf 


then  Kg  becomes 

K F  =  lan  TH 

Thus,  for  high  Q  filters  (Q  >  10)  this  can  be  approximated  by 


U-27 

11-28 


II-29 


11-30 


11-31 


kb  T3 


11-32 


As  expected,  the  higher  the  Q  of  the  filter,  the  narrower  the  bandwidth.  If,  for 
example,  k  is  chosen  to  be  equal  to  four  --  signifying  that  the  sampling  rate  is 
four  times  the  desired  carrier  frequency  and  Q  «  tt  x  100/2  (for  example),  then 
the  digitized  R-C  of  II-28  becomes 
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H(Z)  = 


1  +  Z 


1.005  +  . 995  Z 


rz- 


A  plot  of  this  filter  within  the  rejection  band  is  shown  in  Figure  II-2. 

Note,  once  again  that  this  filter  has  been  plotted  as  a  function  of  angle. 
Since  it  was  assumed  that  i  IT  =4,  then  w  T  =  tr/4  or  90°.  Thus  the  normalized 

S  O  c 

center  frequency  of  this  bandpass  filter  is  90°. 

The  digital  LP-BP  and  LP-BS  transformations  given  in  11-21  and  11-22 
have  interesting  properties  when  compared  with  the  LP-LP  transformation  of 


S 


<1  -  Z'1) 
K-  (1  +  Z*1) 

GJ 

u 


11-34 


obtained  by  combining  the  bilinear  transformation  with  S  — *■  S/uj^. 

Consider  11-21  when  the  ratio  of  sampling  rate  to  carrier  frequency 

f  /r  is  4:1.  In  this  case  cosiU  T  =  0  and  11-21  becomes 
s'  o  o 


1  <1  +  Z"2) 

kE‘  (i  -  z"z) 


Similarly  the  LP-BS  transformation  becomes 


H-35 


_  y  (1  -  Z'2) 


11-36 


Comparing  11-34  and  11-35,  one  notes  that  if  K—  is  replaced  by  Kw-  and  Z  is 

U>  D 


replaced  by  -Z  then  11-34  ;s  converted  to  11-35.  Thus,  one  has  a  procedure 
to  convert  a  lowpass  digital  filter  to  a  bandpass  digital  filter  without  changing 
the  magnitude  of  any  of  the  coefficients.  The  constraints  are  that  the  upper 
frequency  cutoff  of  the  lowpass  filter  is  made  to  coincide  with  the  digital  band¬ 
width  and  the  4:1  ratio  be  maintained.  Thus  a  digital  lowpass  filter  of  the  form 
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11-37 


i 


N 

y  a  Z-n 


«LP<Z»  -  T 


Z  \Z‘" 

n  =  0 


is  converted  to  a  digital  bandpass  filter  of  the  form 

N 

,n  „-2n 


Z  bn(-l)n  Z 

•— i  n 

n  =  0 


n  ~-2n 


11-38 


Similarly,  the  LP-BS  transformation  is  obtained  by  replacing  Z  by 

Z~2  and  K—  by  l/Kj£  •  As  an  example  of  the  latter,  consider  applying  11-34 
u 

to  the  analog  R-C  filter  of  11-24.  Making  the  substitution  yields 


h,p<z>  ■  — 


1  +  Z 


-1 


-1, 


,-l. 


(■j ^T-)d  -  Z  )  +  (1  +  Z  ) 

u 


11-39 


Making  the  substitution  noted  above  yields 


HBS<Z>  = 


1  +  z 


-2 


K-g  ( 1  -  Z'Z)  +  {1  +  Z*2) 


11-40 


hbs<z> 


1  +  z 


( 1  +  K  fr  )  +  (1  -  Kwr-)Z 


11  ■  u 


Equation  11-41  reduces  to  the  notch  of  11-33  when  K  .  005  which  is  in  the  ex¬ 
ample  . 

The  above  procedures  illustrate  the  use  of  a  digitized  LP-BP  and  BS 
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transformations .  As  noted  above,  when  the  ratio  of  sampling  to  carrier  fre¬ 
quency  is  4:1,  bandpass  filters  can  be  obtained  without  changing  the  magnitude 
of  the  coefficients.  It  is  to  be  noted  that  the  highest  negative  power  of  Z  is  twice 
what  it  was  in  the  lowpass  case,  implying  twice  the  number  of  delay  elements. 


C.  BANDPASS  DIGITAL  FILTERS  -  A  SHIFTING  TECHNIQUE 

The  bandpass  design  procedure  discussed  in  the  foregoing  utilized  the 
LP-BP  transformation  in  conjunction  with  the  bilinear  transformation.  The 
analog  LP-BP  transformation  is  of  a  form  which  yields  a  realizable  (in  terms 
of  R's,  L's  and  C's)  analog  filter.  One  need  not  be  constrained  in  this  manner 
in  the  design  of  digital  filters.  A  technique  will  now  be  described  which  parallels 
that  of  a  heterodyning  procedure  in  analog  signal  theory. 

If  one  has  the  Fourier  transform  pair 

f(t)  *-*  F(w)  11-42 

Then 

F(w  -  u>  )  +  F(w  +  w  ) 

g<t)  -  f(t)  cos  wot  — - 2__ - II-43 


Thus  il  f(t)  were  a  lowpass  signal,  then  g(t)  is  a  bandpass  signal  centered  at 
Consider  the  application  of  the  above  to  the  lowpass  digital  filter 

N 

V  a  Z_n 


hlp(Z)  •  ^zr- 


11-44 


V  b  Z’n 
n 

n  -  0 


By  dividing  out  the  denominator- ,  11-44  can  be  written  as 


oo 

HLP(Z>-  I  hnZ'n 

n  =  0 


n-45 


s 
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Thus,  the  numerator  coefficients  of  the  BP  filter  are  a  function  of  both  the  nu¬ 
merator  and  denominator  coefficients  of  the  lowpass  filter.  It  should  also  be 
noted  that  the  bandpass  coefficients  are  once  again  a  function  of  the  normalized 
carrier  frequency  and  that  11-48  and  11-49  are  in  the  form  of  a  convolution  sum 
suggesting  a  synthesis  procedure,  which  will  be  discussed  in  the  following 
sections. 

Perhaps  the  greatest  application  for  this  technique  lies  in  the  design  of 
symmetrical  bandpass  filters  --an  application  which  will  b  discussed  at  length 
in  the  next  section.  The  above  application  is  due  to  the  fact  that  the  design  pro¬ 
cedure  was  accomplished  through  a  translation  of  the  frequency  characteristic 

to  the  left  and  right  by  w  .  The  LP-BP  transformation  discussed  previously 

2 

had  properties  of  geometric  symmetry 

As  in  the  case  of  the  LP-BP  transformation,  a  simplification  is  obtained 
for  a  4:L  ratio  of  sampling  frequency  to  center  frequency.  For  this  case,  11-50 
and  11-51  become 

M 

A  =  (- l)n  V  {_i)r  b  a  11-52 

n  ^  r  n  -  r 

r  =  0  n  even 

M 

B  =  (-l)n  V  (-  l)r  b  b  11-53 

n  '  '  — /  r  n-r 

r  =  0 

A  =  B  =  0  n  odd 

n  n 

Illustrative  examples  comparing  the  above  and  the  LP-BP  transformation  will 
be  discussed  in  the  next  section.  The  two  approaches  presented  provide  for 
convenient  synthesis  of  bandpass  (and  bandstop)  digital  filters  in  terms  of  their 
critical  parameters. 
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D.  BANDSTOP  DIGITAL  FILTERS  -  A  FEEDBACK  TECHNIQUE 

Consider  a  bandpass  digital  filter  given  as  HBp(Z).  Then  form  a  new 
digital  filte-  given  by 


hbs(z> 


1 

i  +  K  hbp(Z) 


11-54 


where  K  is  a  constant  adjusted  to  be  much  greater  than  the  maximum  gain  of 
HfipfZ).  Now,  assuming  that  the  gain  of  .iBp(Z)  is  unity  in  the  pas-band  then 
Hbs(z)  is  given  by 


Hns(Z)  1  +  K 


K 


fin  the  pass>-“nd' 
1  ofHBp(Z) 


11-55 


In  the  region  where  HBp{Z)  is  small  (KHBp(Z)  «  1),  Hgg  is 


HBS<Z>  «  1 


(In  stopband! 

of  hbp(Z)  J 


11-^6 


Thus,  the  alternation  in  the  stopband  is  controlled  by  K  and  the  passband  gain  is 
unity. 

Consider  the  digital  filter  coefficients  generated  by  this  technique.  Let 


HBp(Z) 


N(Z) 

Tn?) 


11-57 


Then  11-54  becomes 


hbs<z> 


D(Z) 

wrofiy 


11-58 


Thus,  the  numerator  and  denominator  coefficients  are  given  respectively  by 


a  =  b 
n  n 


11-59 


b'  =  b  +  Ka 
n  n  n 

The  primed  coefficients  are  those  of  the  resulting  bandstop  filter.  The  unprim¬ 
ed  coefficients  are  those  of  the  bandpass  filter  .  If  this  technique  is  combined 
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with  the  lowpass  to  bandpass  technique  discussed  in  section  B,  then  11-59  be¬ 
comes 


a'  =  (-l)n  b 
2n  n 

bL  =  (-l)n  b  +  K( - 1  )n  a 
2n  n  n 

a1  =  b'  =0;  q  odd 

q  q 


n-6o 


The  unprimed  coefficients  are  now  those  of  the  lowpass  filter. 

As  an  example  of  this  technique,  consider  a  normalized  fourth  order  low- 
pass  Chebycheff  filter  with  a  l/2  db  ripple.  The  coefficients  for  this  analog 
filter  are  tabulated  in  many  analog  filter  design  handbooks  and  are  given  as 

TABLE  VI 

Analog  Fourth  Order  Chebycheff  LP  Filter 
( 1/2  db  ripple) 


Ao=  1 

Bo  = 

.  379 

A,  =  0 

Bi  = 

1.025 

A2  -  ® 

B2 

1,717 

> 

II 

o 

B3 

1.  197 

> 

II 

o 

B4 

1 

The  coefficients  for  the  lowpass  digital  filter  with  a  20°  cutoff  frequency  (18:1 
ratio  of  sampling  rate  to  cutoff  frequency)  become  (using  the  expansion  proced¬ 
ure  associated  with  Table  I) 
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TABLE  VII 


Digital  Lowpass  Chebycheff  Filter 
(20°  cutoff) 


aQ  =  9. 667  x  10 

at  =  3.867  x  10‘3 

a2  =  5. 800  x  10*3 

a3  =  3. 867  x  10'3 

a  .  =  0.  667  x  10'4 
4 


bi  =  -4.410 

b2  =  5.896 

b3  =  -3.588 

b  .  =  .837 

4 


Ut.lizing  11-60,  the  bandstop  filter  coefficients  are  generated  with  K  chosen  as 

4 

K  =  10  (an  80  db  attenuation).  These  digital  filter  coefficients  are  then  given  as 


TABLE  VIII 

Bandstop  Chebycheff  Filter 
(carrier  frequency  =  90°) 


Bandwidth 

0 

o 

CM 

11 

ao  = 

1. 271 

bo  = 

10. 

938 

11 

4.  410 

b2 

34. 

260 

a4 

5.896 

ii 

-  TJ* 

J3 

63. 

899 

2.  588 

b6 

35. 

082 

a8 

.  837 

b8 

10. 

504 

An  interesting  feature  of  this  technique  is  that  the  resulting  "a"  and  "b"  coeffi¬ 
cients  (Table  VIII)  are  approximately  the  same  order  of  magnitude,  although  the 

-4 

lowpass  filter's  numerator  coefficients  (Table  VII)  arc  approximately  10  of  the 


denominator  coefficients.  In  fact,  it  can  be  shown  that  as  the  lowpass  filter's 
bandwidth  is  made  smaller,  the  "a"  coefficients  also  get  smaller.  Thus,  one 


can  design  a  notch  or  band  rejection  filter  with  increasing  notch  depth  as  the 
notch  width  gets  narrower  and  still  maintain  a  small  dynamic  range  for  the 
filter  coefficients.  The  bandpass  and  bandstop  filter  coefficients  of  Tables  VII 
and  VIII  are  plotted  in  Figures  II-3  and  II-4. 

The  technique  utilized  here  is  analogous  to  placing  the  bandpass  filter  in 
the  feedback  portion  of  a  loop  which  has  a  forward  loop  gain  of  K.  If  a  bandstop 
filter  were  placed  in  the  feedback  loop,  the  resulting  filter  would  be  a  bandpass 
filter . 

Various  transformation  techniques  have  been  discussed  for  use  in  the 
generalized  design  of  lowpass,  bandpass  and  bandstop  filters.  The  concept  of 
normalization,  common  to  all  of  the  above,  was  noted  to  be  of  extreme  impor¬ 
tance  in  the  design  of  digital  filters.  Techniques  satisfying  particular  constraints 
as  well  as  various  error  considerations  will  be  discussed  in  the  following 
sections . 
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f  'GuRe  n-3 
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III.  ARITHMETICALLY  SYMMETRICAL  BANDPASS  FILTERS 

A  design  technique  was  investigated  which  yielded  bandpass  digital  filters 
with  equal  ripple  amplitude  characteristics  in  the  passbaad  and  stopband,  having 
arithmetically  symmetrical  cutoff  characteristics.  Filters  having  arithmetic 
symmetry  have  application  in  many  AM  and  FM  problems.  Conventional  analog 
filters  which  have  been  transformed  to  bandpass  have  geometric  symmetry  prop¬ 
erties,  The  approach  taken  was  to  design  an  equiripple  or  elliptic  analog  filter 
at  lowpass  and  then  convert  this  filter  to  a  lowpass  digital  filter  through  the  use 
of  the  bilinear  transformation.  Finally,  these  filters  were  converted  to  a  band¬ 
pass  digital  filter  while  preserving  the  desired  properties.  Errors  and  constrain¬ 
ing  relations  among  the  parameters  are  dis<  ussed.  The  foregoing  procedures 
will  now  be  detailed. 

A.  DESIGN  PROCEDURE  FOR  ELLIPTIC  ANALOG  FILTERS 

It  was  desired  to  formalize  a  design  procedure  to  be  used  for  elliptic 
analog  filters.  As  discussed  by  Calahan[^  this  procedure  can  be  formalized  as 
follows : 


Consider  the  lowpass  characteristic  as  shown  below 

1 

!  K 


Figure  III  -  1 
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The  equinpple  characteristic  is  such  that 


a.  at  w  =  0  T(u)  =  K 


b .  at  w  =  1  |  T  (o) 


1 2  K 


c  .  at  w  =  u  |  T(cj) 


2  _  K 


(u  >  w  defines  the  stopband;  w  <_  1  defines  the  passband.  ) 

With  the  above  definitions,  the  following  procedure  is  used 

2  2 

1.  Calculate  m  =  l/oj“  nn  =  (t,/en) 

'  c  r  2 

2.  Determine  the  order  of  the  filter,  N,  by  choosing  N  as  the  smallest 


integer  satisfying 


Kfm'j)  K(m) 
K(mJ)  K(m') 


where  =  !  -m,  m'^  =  1  -rn'  and 


K(m)  =  j 


0  (1  -m  sin  x) 


which  is  the  real  quarter  period  of  the  elliptic  integral  of  the  first 


N  can  be  approximated  for  values  of  m  near  1  and  m'  near  zero  by 
N  >  < 2/nZ )  in(4e2/tl)  in  8/(«c-  1) 

3.  Determine  the  zeros  and  poles  of  the  transfer  function  from  the 


formulas 


2  -  _  _ ~  J _ 

T  m sn(£-  K(m),  m) 


J  sn(  t  —  K(m)  +  j  r„ ;  m) 


where  r  =  2,  4,  •  •  •  ,  N  -  1  N  odd 
r  =  l,3,«**,N-l  N  even 


r 

o 


K(m)  bc*1  {■—  jm'  ) 
_ £  1  1 

NK(m') 


sinh'1  <-L) 
K(m')  £1 


for  m'  near  zero. 

The  term  sn(-)  is  the  elliptic  3ine  and  the  term  sc"*(4>:m)  is  the 
elliptic  arc  tangent. 

As  an  example  of  this  procedure,  consider  the  following  filter  specifica¬ 
tions 


Passband  -  2db  ripple  0  <  u  <  1 


Stopband  -  25 db  min.  1.33  <  w  <  oo 
For  this  example,  the  above  parameters  become 
rn  =  .  566  K(m)  =  1. 909 


Cj  =  .765 


v  =  .438 

o 


e  2  =  161 


N  =  2.  98  (choose  N  =  3) 


*  stjl.47  Pj  =  -.463  p2  =  129  +  j.  962 


With  the  above  pole  zero  configuration,  the  resulting  transfer  function  becomes 


C(S)  = 


S  +  2.  16 

- - - 

S  +  .  7 2 1  S  +  1.061S  +  .437 


m-i 


Using  the  above  approach,  several  lowpass  elliptic  filters  will  be  synthesized 
and  then  digitized  using  the  bilinear  transformation. 
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B.  DIGITIZED  LOW  PASS  ELLIPTIC  FILTERS 

The  general  pole -zero  configuration  for  a  lowpass  elliptic  filter  is  shown 

below. 


The  pattern  shown  is  for  a  seven  pole,  or  seventh  order  elliptic  filter.  The  pole 
positions  for  these  filters  lie  on  an  elliptical  contour.  The  values  for  these 
poles  and  zeros  were  determined  so  ifiat  they  would  satisfy  the  following 
specifications : 

Passband  Spec .  <  ,  5  db  r  ipple  0  <  w  1 

Stopband  Spec  .  60  db  min.  1 . 25  <  ^  co 

A  seventh  order  filter  with  the  following  characteristic  was  found  to 
satisfy  these  requirements. 
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TABLE  I 


p0  =  . 3889  z2  =  t  j  2. 391 

Pj  =  .0377  t  j  1.0120  z  =tjl.2711 

Pj  =  .  1399  t  j  .8939  z^r+j  1.4709 

p5  =  . 2932  +  j  .573  1 

Applying  the  bilinear  transformation  to  this  filter  with  a  samplinp  rate 

equal  to  4  times  the  cutoff  frequency  -  -  or  a  normalized  cutoff  at  w  T  =  90° 

(K—  =  tan(45°)  =  1)  yields 
u 


where 


TABLE  II 


aQ  =  . 556  bQ =  12. 252 

aj  =  2. 007  b  =  -14.  381 

a^  =  4. 254  b2  =  30. 397 

a„  =  5.  976  b.,  =  -26.  018 

3  j 

a  .  =  5.  976  b.  =  24.  095 

4  4 

ae  =  4.  254  b„  =  -13.  421 

a ,  -  2.007  b,  =  5.  917 

o  o 

.•7  =  .  556  l»?  =  -  1,519 
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A  plot  of  this  filter  is  shown  in  Figure  III  -2. 

A  second  design  was  carried  out  for  a  filter  with  greater  attenuation  in 
the  stopband  but  a  larger  transition  region.  The  specifications  on  this  filter 
were 

Passband  Spec.  <  .  5  db  ripple  0  <  u  <  1 
Stopband  Spec.  90  db  min.  2  <  to  <  1 

A  seventh  order  filter  with  the  following  pole  zero  pattern  was  found  to 
satisfy  these  requirements. 

TABLE  III 

P0  =  -32o  z2  =  t  3  4-  354 

p  =  .  057  +  j  1. 016  z4  =  t  j 2. 044 

p^  =  .  17.'  t  j  .841  *6  =  t  j  2.  490 

Pj.  =  .  278  t  j  .  486 

Applying  the  bilinear  transformation  to  this  filter  with  a  normalized  cutoff  at 

«  T  =  90°  yields 
u  J 

7 

V  a  Z'n 

— ■  n 

H(Z)  =  -  III -3 

V  b  Z’n 

— -  n 

n  -  0 

where  the  coefficients  a  and  b  are  as  given  in  Table  IV. 

n  n  ** 


TABLE  IV 


aQ  =  . 744  bQ  =  11. 077 

a:  =  4.  075  bj  =  -16.  234 

a2  =  10.  465  b2  =  28.  883 

a3  =  16.  172  b3  =  -27. 885 

a  =  16.  172  b  =  23. 359 

4  4 

a5  =  10.  465  b5  =  -13.  562 

a6  =  4.  075  b6  =  5.  630 

a_  =  .  744  b_  =  -  1.  37 1 

<  ( 

A  plot  of  thiB  filter  is  shown  in  Figure  III-3. 

A  third  elliptic  filter  synthesized  was  that  discussed  in  Section  A  with 

the  normalized  cutoff  frequency  at  =  20°.  The  characteristic  for  this 

filter  is 


V  b  Z*n 
U  n 

n  =  0 


III  -  4 


where 


TABLE  V 


aQ  =  .  188  bQ  =  1 .  162 

a j  =  -.  14  1  bj  =  -3.  087 

a2  =  - ;  141  b2  =  2.  847 

a3  =  .  188  b3  =  -  . 903 
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F  \  GURE.  HI -4!  LOW  PASS  ELLIPTIC  FILTER 


A  plot  of  this  filter  is  shown  in  Figure  III-4. 

A  technique  which  converts  these  lowpass  digital  filters  to  bandpass 
filters  will  now  be  discussed. 


C.  APPLICATION  OF  THE  SHIFTING  TECHNIQUE 

It  was  initially  decided  to  utilize  the  filter  of  Figure  III-4  for  the  synthe¬ 
sis  of  an  arithmetically  symmetrical  bandpass  filter.  The  technique  used  is  the 
shifting  technique  described  in  Section  II  and  qualitatively  consists  of  shifting 
the  lowpass  characteristic  to  the  left  and  to  the  right  by  an  amount  (in  degrees) 
equal  to  the  normalized  center  frequency.  As  discussed,  this  is  simply  the 
digital  counterpart  of  multiplying  a  time  waveform  by  cos  co^t  and  observing  the 
resulting  spectral  plot. 

The  required  equations  for  this  transformation  are  given  in  Section  II  as 


M 

A  =  V  a  b  cos  w  T  (n  -  2r) 

n  U  n-r  r  o 

r  =  0 
M 

B  =  T  b  b  cos  u>  T  (n  -  2r) 

n  Lt  n-r  r  o 

r  =  0 


III- 5 


where  the  "small"  a's  and  b's  are  the  numerator  and  denominator  coefficients 
of  the  lowpass  filter  and  the  "capital"  A's  and  B’s  are  the  corresponding  band¬ 
pass  coefficients.  From  this  equation  it  can  be  seen  that  the  bandpass  coef¬ 
ficients  are  obtained  by  a  weighted  convolution  (or  correlation)  of  lowpass 
coefficients,  where  the  weighting  function  depends  on  the  desired  carrier  fre¬ 
quency. 

Utilizing  these  equations  with  a  normalized  center  frequency  of 

w  T  =  90°  yields 
o 
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H(Z)  = 


E  a2„z-2" 

n  =  0 

E  b2„  z‘2n 

n  =  0 


m-6 


where 


TABLE  VI 


> 

© 

II 

.  219 

Bos 

1. 351 

A2  " 

.  063 

B2 

2.  910 

> 

n 

.  053 

B4  = 

2.  528 

A6 

-.  170 

B6  = 

.  816 

A  plot  of  this  filter  is  shown  in  Figure  III  -  5 . 

A  detailed  computer  analysis  of  the  magnitude  characteristic  of  this 

filter  shows  that  this  filter  is  arithmetically  Symmetrical  about  90°  with  no 

symmetry  error.  A  detailed  comparison  of  the  lowpass  elliptic  filter  in  Figure 

III  -  4  and  the  bandpass  filter  in  Figure  III  —  5  shows  that  the  passband  error 

between  these  two  filters  is  less  than  1/2%.  These  errors  are  due  to  the  tails 

* 

of  the  shifted  characteristic  extending  back  into  the  passband  of  the  filter. 

Based  on  the  above  results,  it  was  conjectured  that  the  symmetry  error 
can  be  reduced  to  zero.  The  following  derivation  will  show  this  to  be  so  for 
particular  values  of  the  normalized  (with  respect  to  the  sampling  rate)  carrier 
frequency. 


In  comparing  Figures  III  -  -4  and  III  -  5  F  the  factor  of  2  from  Eq.  11-47 
must  be  included. 


45 


.OH  *  If  Of>£  7  nOL  SSVd  QMV9 

.O^  -  =  331^33 

»JL<?18,+  »-g-PZfi-Z  +  ,-EO>6'  2  +  l^ei  ^ 

9.iOtl--».ltSO'+,.J4<50-+  b\Z'  ~~  '  ' 


l &URE  UI.-5:&/kMD  PA.SS  ELLIPTIC  FtLT 


D.  SYMMETRY  ERROR  AND  CHOICE  OF  CENTER  FREQUENCY 

Consider  the  lowpass  filter  to  be  expressed  as  a  function  of  angle,  that 
is  FLp(6),  where  9  =  coT  and  w  is  the  digital  frequency  and  l/T  =  fg,  the  sam¬ 
pling  rate,  Now  due  to  the  periodicity  of  digital  filters  one  has 

FLp(6)  =  FLp(6  +  2t)  IF-  7 

Because  of  symmetry  requirements 

I  FLP<0)  I  =  lFLP(_e)1  m-8 

Therefore 

|FLp(-Tr  +  6)|  =  |FLp(rr  -  9)  |  III- 9 

Now 

FlP(-tt  +  0)  =  FLp(-tr  +  8  +  2  it)  =  FLp(w  +  9)  III  - 1 0 

or 

|FLp(u  +  0)|  =  |  F Lp(Tr  -6)|  HI- 11 

From  the  discussion  in  Section  II,  it  can  be  shown  that  the  bandpass  filter  result¬ 
ing  from  utilizing  this  shifting  technique  yields  a  bandpass  filter  Fgp(9)  given 
as  (neglecting  a  scale  factor  of  2); 

fBP|a>  *  FLP(0  +  V  +  FLP(6  ’  V  IU-12 

where  6  =  oj  T  and  w  is  the  desired  center  freouencv. 

o  o  o  * 

Now  consider  the  response  of  the  bandpass  filter  at  a  distance  (in  angle) 
away  from  0  .  Thus, 

FBP(9o  +  &  =  FLP(29o  +  ^  +  FLP(4>) 

in-13 

FBp<eo‘*>sFLP<2eo-+>  +  FLP(-+) 

In  order  for  symmetry  to  exist  at  bandpass,  Eqs.  Ill  - 1  3  would  have  to  be  such 


that 


i  FbP(9o  +  ^  I  =  I  FBP(6o  ‘  ^  I  m‘14 

Moreover,  :n  order  for  there  to  be  no  shifting  error  from  lowpass  to  bandpass, 
one  requires  that 


FBP(0O  +  *)  -  flp(*) 

fbp(0o  -  ♦)  -  rLP<-«.) 

(or,  the  above  modified  by  a  scale  factor). 

As  can  be  seen  by  comparing  Eqs.  Ill  - 1 5  and  IH-13, 


III- 15 


this  latter  error  is 


LP-BP  SYM.  ERROR  =  F.  D(Z0  +  <J>)  HI- 16 

J _j  i  o 

or 

FLP«Zflo*« 

Consider  the  case  where  0q  =  ir/Z.  This  choice  of  center  frequency  is  equiva¬ 
lent  to  a  sampling  rate  of  four  times  the  center  frequency.  That  is 


0  =  u>  T  - 

o  o 


2  it  T 

o 

~1 


s  s 

For  the  above  choice  of  0  ,  Eqs.  Ill- 13  become 

FBP(6o  *  *>  =  flp<"  +  *>  +  flp(*> 

FBP<6o  ‘  41  =  FLP^n  -  41  ♦  FLP(‘‘! 
From  Eqs.  Ill -9  and  III  - 1  1,  Eq.  Ill- 18  becomes 


III- 17 


III- 18 


FBp(eo  +  4>)  =  Mj  <  Aj  +  M2  <  A2  =  K  <  Q  III- 19 

FBP^9o  '  <>)  =  Mi  <  +  M2  <  -&2  ~  FBP*So  +  ^ 

Thus  for  a  4  to  1  ratio  of  sampling  rate  to  carrier  frequency  there  is  no  sym¬ 
metry  error  at  bandpass.  The  error  between  the  lowpass  filter  and  the 
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bandpass  filter  is  given  by 


where  the  filter  coefficients  become 


TABLE  VII 


A0  =  .2  1875370  B0  =  1.35144966 

Aj  =  .74456616  Bj  =  7.17729147 

A2  =  .80670596  B2  =  16. 1488383 

A3  -  08252552  B3  =  19. 6781 1 17 

A .  =  -.85454586  B .  =  13.6836573 

4  4 

A,  =  .66294855  B5  =  5.14441316 

A,  17000589  B,  =  .8162383720 

b  o 

A  plot  of  this  bandpass  filter  is  shown  in  Figure  III-6.  A  comparison  of  the 
computer  print-outs  for  this  filter  and  the  equivalent  lowpass  filter  bears  out 
the  foregoing  analytical  conclusions. 

It  is  of  interest  to  note  the  recursion  equations  which  result  from  the 
above  choices  for  0^.  In  general  these  recursion  equations  are 

M 

A  =  ^  a  b  cosfo  (n  -  2r)l 

n  t-  n  -  r  r  o  J 

r  =  0 


B 

n 


M 

^  b  b  cos  f©  (n-2r)1 
—  r  n-r  (.  o'  '} 

T  -  0 


III-24 


where  G  =  u>  T. 
o  o 

For  0  =  it  (a  two-to-one  ratio),  these  equations  become 

o 
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FI  eueEHL-^:&AWD  PASS  EL  L  IP  TIC  FILTER 


JIT-25 


A 

r. 


M 


(-1) 


n 


V 


a  b 
n-r  r 


r  =  0 


B 

n 


(-1) 


n 


M 

V 


r  -  0 


b 

n  -  r 


Thus,  the  bandpass  filter  coefficients  (A^,  B^)  are  the  result  of  convolving  the 
appropriate  lowpass  filter  coefficients  (a^,  b^)  with  each  other.  Since  a  digital 
filter  performs  the  operation  of  convolution,  this  suggests  the  following  auto¬ 
mated  procedure  to  generate  the  bandpass  coefficients  once  the  lowpass  filter 
has  been  implemented.  If  one  assumes  that  the  lowpass  filter  has  been  synthe- 
sized  using  the  direct  form  and  the  "x  bank1'  or  feed  forward  section  is 
impulsed,  then  the  outputs  are  the  "a"  coefficients.  If  these  outputs  are  fed 
into  the  "y  bank"  or  feedback  section,  its  outputs  are  the  "A"  or  numerator 
coefficients  of  the  bandpass  filter.  The  "B"  coefficients  can  be  obtained  in  a 
similar  manner  using  only  the  feedback  bank.  This  procedure  appears  to  be  of 
advantage  in  that  the  lowpass  coefficients  can  be  altered  and  the  lowpass  fiber 
itself  used  to  generate  the  bandpass  equivalent. 

The  equations  analogous  to  I  IT  -  5  for  0Q  =  tt  /  2  (a  four-to-one  ratio)  are 


A  =  (cos  n  -J—  )  ^  { -  1  )r  a  b 

n  2  -j  n-r  r 

r  -  0 
M 

B  =  (cos  n  -Z-)  V  ( _  l  )r  b  b 

n  2.  '  n-r  r 

r  -  0 


III  -  26 


This  approach  is  not  limited  to  any  one  implementation. 
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A  similar  implementation  i.  an  be  utilized  for  the  above.  It  should  also 
be  noted  that  and  B^  are  both  zero  lor  n  odd.  Thus,  although  the  total  delay 
required  for  both  bandpass  filters  are  the  same,  only  half  the  multiply -add 
operations  are  required  when  a  four-to-one  ratio  (f g/ ^  =  4}  is  used, 

E.  APPLICATION  OF  THE  LP-BP  TRANSFORMATION 


In  Section  II,  the  lowpass  tc  bandpass  transformation 

1  -  2(cos  w  T)  Z'1  +  Z‘Z 

S  _ - 2 - - -  III-27 

K  b  ( 1  '  Z  **> 

was  utilized  in  the  design  of  bandpass  digital  filters.  Mention  wau  also  made 
of  the  simplification  resulting  from  adjusting  the  sampling  ' :  to  be  four  times 
the  carrier  frequency.  This  choice  causes  so  1  -  rr/2  and  Eq .  Ill  -  27  becomes 


1  +  Z 

kb°  - 


-2 


Z 


HI -28 


It  is  desired  to  investigate  the  application  of  the  LP-BP  transformation  to  the 
design  of  arithmetically  symmetrical  filters  for  tne  above  choice  of  normalized 
digital  center  frequency.  This  transformation  is  recognized  as  being  that  of  the 

.  j  -2  -2 

lowpass  bilinear  transformation  with  Z  replaced  by  -  Z  .  Since  -Z  can  be 
written  as 


,  4-  ■  •  2  j ( uj  T  t  -5-  ) 

~-2  ;.J1T  j2uT  J  2  1  tii  in 

-Z  =eJ  e^  =e  III-29 

This  transformation  can  be  viewed  as  a  shift  of  90°  accompanied  by  a  scale 
change  by  a  factor  of  2.  Thus,  with  these  parameters,  a  lowpass  filter 
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i 


nLP(/)  "  ~TT 


N 

V  a  Z'n 
n 

n  =  0 


111-30 


7  b  Z-n 
n 

n=  0 


become  s 


N 


HBP(Z)  ‘  N 


l  a(- DnZ 
n  =  0 


n  ,  -2n 


III  -  3  1 


n  _-2n 


l  b„'-»"z 

n  -  0 

Now,  consider  the  bandpass  filter  evaluated  at  (0  +  tt/2)  =  uT.  Equation  IU-31 
become  s 

N 


HBP(Z)  '  N 


7  a  (-l)n  c-i*  n  e’2jn0 
Lj  n 

n  =  0 


IU-32 


T  b  <-l)n  e'jlT  n  e’j2n9 
<  n 

n=  0 


N 


y  _  9 


a  e 

n 


hbp(Z>  =  ^r 


III- 3  3 


b  e 
n 


■2  j  n  0 


n  --  0 


This  is  recognized  as  being  identical  to  the  lowpass  filter  evaluated  at  2  0. 
Similarly,  Hgp(ir/2  -  0}  =  H^p(-20).  Thus,  this  technique  yields  an  arith¬ 
metically  symmetrical  lowpass  filter  with  no  symmetry  error  either  at  band¬ 
pass  or  in  the  shift  from  lowpass  to  bandpass.  Two  examples  of  this  technique 
are  a  4th  order  Chebycheff  lowpass  filter  and  a  10th  order  Butter  worth  lowpass 
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i 

i 

1 

« 

! 


\ 

! 

] 


1 


filter  that  have  been  shifted  to  bandpass.  '1  he  coefficients  for  these  filters  are 


TABLE  VIII 
10th  Order  iiutterworth 


Lowpa 

ss  Coefficients 

Bandpas  s 

Coefficients 

tan 

(?0T/Z)  = 

1 

w  T 
o 

=  tr/2 

Numerator 

Denominator 

Numerator 

Denominator 

ao 

i 

b0 

345.  252 

A0  " 

i 

B-  * 

KJ 

345. 252 

31  = 

10 

bl  = 

-4  x  10~J 

A,  - 

far 

-  10 

BZ  = 

4  x  10* 

a2  = 

45 

b2 

462.771 

A4 

45 

B4- 

462.771 

a3  = 

120 

b3 

8  x  10'3 

A6 

-120 

B6 

-8x10 

3  4  = 

210 

b4 

18b. 270 

A8  = 

210 

B8  = 

188. 270 

a5  = 

252 

b5 

-8  x  10"5 

A10  " 

-252 

L10  = 

8  x  10* 

a6  = 

210 

b6 

26. 613 

A1  2  = 

210 

B 1  2  = 

26.613 

a.,  - 
( 

120 

b,  = 

1 

-9  x  10*3 

A  1 4  = 

-120 

B 1 4  ~ 

9  x  10* 

a8 

45 

b8 

1.  093 

Al6 

45 

B16 

1. 093 

3  9  = 

10 

b9 

-4  x  10'3 

A 18  * 

-  10 

B 1 8  : 

=  4x10' 

a10  = 

1 

b10  = 

.  024 

A20  = 

1 

B20  = 

.  024 

(A  =  0; 
n 

n  odd) 

<Bn  = 

0;  n  odd) 
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T AI1LL  IX 


4th  Order  Chebychcfl 

Lowpass  Coefficients  Bandpass  Coefficients 


tan(wT/2)  =  , 

.  17633  (20° 

cutoff) 

^  T  = 
o 

tr/2 

Numerator 

Denominator 

Numeral  or 

Denominator 

ao  =  1 

bo- 

5.319 

Ao  = 

i 

B0 

=  5.319 

a,  =  4 

b!  = 

-2. 828 

A2 

4 

B2 

=  2.828 

a2  ”  6 

b2 

4.  841 

A4 

6 

B4 

*  4. 140 

a3  =  4 

b3 

-2.  140 

A6 

4 

B6 

=  2. 140 

II 

b4  = 

.  870 

A8  “ 

1 

B8 

=  .  870 

(An  =  0;  n  odd)  (Bn  =  0;  n  odd) 

These  filters  are  plotted  in  Figures  III -7  and  III  -  8 . 

It  is  of  interest  to  note  that  the  scale  change  caused  by  this  transfornr.a  - 

tion  yields  a  bandwidth  at  bandpass  that  is  i  ne  half  that  of  the  lowpass  filter. 

Thus,  for  the  4th  order  Chebycheff  filters  the  bandwidth  is  +  20°  at  lowpass  and 

+  10°  at  bandpass.  Of  course,  since  these  bandwidths  are  normalized  to  the 

sampling  rate,  the  actual  bandwidth  can  be  conti  oiled  by  this  rate.  It  nould 

also  be  noted  that  the  filter  coefficients  generated  are  obtained  directly  from 

the  lowpass  design.  By  contrast,  the  shifting  techr.ique  used  in  the  foregoing 

yields  a  more  complicated  structure.  For  example,  the  third  order  elliptic 

filter  of  Table  V  had  bandpass  coefficients  using  the  shifting  technique  as  gi'  en 

.  - 1  -2 

below.  The  coefficients  obtained  through  the  substitution  Z  -  -Z  are  shown 
in  parenthesis  in  Table  X. 
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lOT  ORDER  BANDPASS  BOTTER WOR.TM 


TABLE  X 


Elliptic  Filter  Coefficients 


Ao  = 

.219  (.  188) 

Bo  = 

1.351  (I.  162) 

A2  ~ 

.  063  (.  141) 

B2  = 

2.  910  (3. 087) 

A4 

.  053  (-.  141) 

B4 

2.  528  (2.  847) 

A6 

-.  170  (.  188) 

B6  ~ 

.816  (  .  903) 

Although  the  coefficients  are  similar,  those  of  the  shifting  technique  are  gener¬ 
ated  through  a  convolution  formula,  whereas  those  of  the  present  technique,  by 
a  simple  sign  change.  Thus,  if  it  is  possible  to  use  a  4:1  ratio,  the  LP-BP 
procedure  is  preferable.  For  other  ratios,  however,  the  LP-BP  procedure 
does  not  yield  satisfactory  results. 

Two  approaches  to  the  synthesis  of  bandpass  digital  filters  having  arith¬ 
metic  symmetry  properties  have  been  discussed.  The  shifting  technique  was 
found  to  be  of  general  applicability  and  was  shown  to  yield  negligible  errors 
under  appropriate  choices  of  the  ratio  of  sampling  rate  to  desired  carrier 
frequency.  The  equations  relating  the  coefficients  of  the  bandpass  digital  filter 
to  those  of  the  low  pass  digital  filter  were  shown  to  be  in  a  form  convenient  for 
"on  line"  synthesis.  Lastly,  the  LP-BP  transformation  was  shown  to  be  an 
excellent  (and  simple)  technique  if  the  problem  constraints  allowed  the  use  of 
a  4:1  ratio  of  sampling  rate  to  carrier  frequency. 
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IV.  NONRECURSIVE  DIGITAL  FILTERS 


The  design  ofa  nonrecursive  (i.e.,  finite  memory)  digital  filter  consists  of 
choosing  coefficients  so  that  a  frequency  response  characteristic  of  the  form 

N  _ 

H{ Z)  =  Yj  hn  Z'n  :  Z~l  =  e"jwT  IV-1** 

n  =  0 

.  -1  -jcoT 

has  the  desired  shape.  From  this  equation  it  is  3een  that  with  Z  =  e 

H(Z)  is  represented  as  a  finite  Fourier  series  with  coefficients  hn.  The  overall 

*Since  the  input-output  equation  of  a  nonrecursive  filter  is  given  by 

N 

y  =  y  h,  X  , 

7n  Li  k  n-k 

k=  0 

it  can  be  seen  that  the  output  samples  are  a  function  of  the  present  input  as  well 
as  past  inputs  up  to  the  upper  limit  of  the  summation.  Thus,  N  is  often  refer¬ 
red  to  as  the  memory  of  the  filter. 

**The  link  between  nonrecursive  digital  filters  and  analog  tapped  delay 
line  or  so-called  transversal  fillers  can  be  seen  by  noting  that  the  input-output 
equation  corresponding  to  IV-1  is  given  by 

N 

v  =  y  h,  x  , 

yn  Li  k  n-k 
k=  0 

whereas  that  of  a  transversal  filter  is  given  by 

N 

y(t)  =  Tj  Ak  x(t-kT) 
k  =  0 

Thus,  if  the  output  is  sampled  at  times  t  =  NT  the  h^  and  A^  are  equal  to  one 
another.  The  Ak  are  the  tap  weightings. 
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question  is,  how  does  one  choose  these  coefficients  in  some  optimum  manner, 

or,  from  what  time  function  do  these  coefficients  come?  If  the  nonrecursive 

filter  is  viewed  as  the  finite  memory  counterpart  of  a  recursive  digital  filter 

which  has  been  synthesized  through  a  transformation  technique  applied  to  an 

analog  filter,  the  hn  coefficients  can  be  related  to  the  impulse  response  of  the 

analog  filter.  If  the  standard  Z  transform  were  applied  to  the  analog  filter, 

then  the  h^  can  be  shown  to  be  related  to  the  first  N  +  1  amplitude  samples  of 

the  impulse  response  of  the  analog  filter.  If  the  bilinear  transform  were  applied 

to  the  analog  filter,  the  h^  coefficients  would  be  related  to  the  first  N  +  1 

coefficients  in  the  Laguerre  expansion  of  its  impulse  response.  Still  another 

approach  is  to  expand  the  desired  frequency  response  in  an  infinite  Fourier 

series  and  then  choose  the  first  N  +  1  terms  or  some  modification  thereof  to  be 

h  coefficients, 
n 

Inherent  in  most  of  the  above  procedures  are  two  general  sources  of 
error  between  the  desired  frequency  response  and  that  finally  obtained.  The 
first  error  is  that  generated  by  the  particular  transform  technique  utilized. 

This  error  can  normally  be  controlled  and  contained  within  tolerable  limits. 
After  obtaining  a  satisfactory  design,  this  set  of  coefficients  is  truncated  to 
achieve  the  finite  memory  indicated  in  IV-1.  This  truncation  procedure  is 
equivalent  to  multiplying  the  time  function  represented  by  the  infinite  set  of 
coefficients  by  a  pulse  of  extent  NT.  This  procedure  produces  an  error  equiv¬ 
alent  to  convolving  the  infinite  memory  response  with  a  sin  x/x  fin  frequency) 
form  whose  central  lobe  width  varies  inversely  with  N,  the  memory  of  the 
filter.  While  this  is  true,  it  appears  desirable  to  have  considerably  more  con¬ 
trol  over  this  type  of  error. 

The  nonrecursive  design  technique  to  be  discussed  in  the  following  sec¬ 
tions  include  design  procedures  from  tabulated  frequency  data,  a  digital 
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y  z 


-k 


k=0 


1  -  Z 


rr 


IV -6 


Equation  IV -5  becomes 


R  -  1 
V 

k  =  0 


H(tk)  e 


*k 

-j  2  TT  l  N 

is  =  y  h 

lj  r 
n  =  0 


-j  2  7T  (i  +  n) 

s 

e  x 


-j  2  it  [t  +  n) 


1  -  e 


R 

TT 


1  -  e 


-j  Z  n  {(  +  n) 


I 

TC 

iv  -? 


If  it  is  desired  to  synthesize  the  response  over  a  half-period,  then  fj  =  0, 
R  =  K/Z  (K  even)  and  IV-7  becomes 
K  , 


-j  Z  n  £ 


H(fk)  e 


k  -  0 


k  N 

z  y  h 

r 

n  =  0 


1  -  e 


-j  it  (i  +  n) 


IV -8 


/K 


for  (f  +  n)  odd  and  k  >  ZN,  zero  for  (f  +  n)  even. 

Although  the  above  equations  can  be  used  to  determine  the  required  h^ , 
it  was  noted  that  a  considerable  simpliiication  comes  about  if  the  desired  digital 
filter  can  be  written  in  the  form 


N 

H(Z)  =  y  h  Z~n 

'  ~i  n 

n  c  -N 

For  this  case,  a  development  paralleling  that  discussed  above  leads  to 


IV  -9 


K  -  1 

V 

«__/ 

k  -  0 


— j  Z  TT  £ 


H(fk)  e 


fk 

7“ 

s  _  V 


N 


-j  Z  it  (f  +  n) 


h  e 
n 


*1 

7" 

s 


n  =  -N 


1  -  e 


-j  Ztt  (£  +  n) 


1  -  e 


-J 


U  +  n j 


iq 


IV-10 


Assuming  that(f  +  r^K  is  not  an  integer,  the  term  in  the  square  brackets  is 
equal  to  zero,  except  for  n  =  -l  ,  In  that  case,  the  form  of  IV  -  1 0  reduces  to 
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K  -  1  -  j  2  Tl  j  t  ~ 


h  =  ^  H(f.  1 

- 1  K  *-J  K- 

0 


;  K  >  2  N 


For  synthesis  over  a  full  period  f  =  -f  /2  and  h4  becomes 


IV- 11 


f 


K  -  1  +  j  2  TT  /  -p 


1  “‘J  17  1  -*•  —  Y  H(f,  )  e  3  '  1 


ht  =  TK  H{~  T_)  e 


r  TT  L,  **'‘k' 
k-  1 


+  TK  H(V  )c 


j  tt  f 


IV -12 

Equation  IV- 12  is  recognized  as  being  equivalent  to  a  trapezoidal  integration  of 

fg/2  j  2  IT  f  -J— 

h  =  y-  J  H(f)  e  S  df  IV- 13 

3  -i  /  2 
8' 

where  H(f)  is  defined  on  K  points  over  f  . 

To  obtain  an  expression  in  terms  of  magnitude  and  phase  characteristics 
one  notes  that 


H(-f)  =H*(f)=  A(f)  e'j  9(f) 
where  A{f)  is  even  and  9,f)  is  odd,  then  becomes 


IV -14 


f./2 


hf  =  j-  j  A(l)  cosfe(f)  -  2w  f  ~  1  df  IV  - 1 5 

s  0  s 


or 


hf  =  TT 


N  -  1 

a(0)  +  y  c°s 

k  =  1 


e(V  _  2~  r‘  r 


T  ( 


f  f 

,  .  n  .  ,  s  ,  . ,  s  , 

-1)  A(-j-)  cos  e(-j-) 

IV- 16 


Thus  ,  relationships  have  been  derived  (IV  -  1  2  and  IV  -  1 6)  which  can  be  utilized 
to  design  nonrecursive  digital  filters  from  tabulated  data. 
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B.  A  DIGITAL  IMPULSE  INVARIANT  TECHNIQUE 

It  is  often  of  interest  to  synthesize  a  nonrecursive  digital  filter  which 
approximates  (in  some  sense)  a  recursive  digital  filter.  Thus,  it  is  assumed 
that  the  recursive  digital  design  satisfies  the  filter  requirements  .  However, 
it  is  desirable  to  eliminate  the  feedback  terms  (b^  coefficients)  and  obtain  a 
nonrecursive  form  which  requires  the  same  amount  of  processing.  Thus,  the 
restriction  might  be  such  that  the  synthesis  is  to  be  carried  out  so  that  the 
total  number  of  coefficients  are  the  same  for  both  recursive  and  nonrecursive 
de  signs . 

Consider  the  transfer  function  of  a  recursive  digital  filter  given  as 

N 

V  a  Z  ~n 
n 

-  IV-|? 

>'  b  z’" 

n 

n  =  0 

By  dividing  the  denominator  into  the  numerator  IV-17  can  be  rewritten  as 


hr(Z)  =  L  c„  z'n 

n  =  0 


IV -18 


where  the  C  are  related  to  the  a  and  b  by 
n  n  n 


'0  b 


(a,  -blCo) 


IV- 19 
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or,  in  general  by 


C 


K  "  TT- 
o 


min  [k.M] 

y 

•-J 

j  =  i 


b. 

j  n 

T“Ck-j 

o 


IV  -20 


where  M  +  1  is  the  number  of  denominator  coefficients  and  the  summation  runs 
from  j  =  1  to  the  smaller  of  the  two  values  k  and  M. 

The  nonrecursive  digital  filter  is  then  of  the  form 


hnr<Z) 


N 

y 

—j 

n  =  0 


h  Z 
n 


IV  -21 


The  question  is  --  how  should  the  h's  be  chosen?  Or,  what  is  the  relation¬ 
ship  between  the  h's  and  the  C's  ?  If  the  h's  are  to  be  chosen  to  yield  a 
minimum  mean  square  approximation  to  H_(Z),  then  the  h^  should  be  equal 
to  the  up  to  n  =  N.  The  justification  for  this  is  that  IV-18  is  in  the  form  of 
an  infinite  Fourier  series  and  IV-21  i3  in  the  form  of  a  truncated  Fourier 
series.  Furthermore,  if  H^j^(Z)  is  to  approximate  a  recursive  filter  synthe¬ 
sis  through  the  use  of  the  bilinear  transformation,  then  the  desired  shape- 
should  first  be  synthesized  with  a  bilinear  transformation  and  then  the  C's 
of  Eq.  IV-21  obtained  through  the  use  ot  Eq .  IV-20.  The  conventional 
approach  is  to  approximate  the  desired  shape  directly  with  a  Fourier  series 
and  then  truncate  after  a  given  number  of  terms  (equal  to  N  +  1  in  the  exam¬ 
ple). 

A  comparison  of  the  various  approaches  is  shown  in  Figure  IV-  1.  The 
ideal  filter  to  be  approximated  is  a  rectangular  filter  whose  bandwidth  (cut¬ 
off  frequency)  is  equal  to  (1/4)  of  the  sampling  rate  (i .  e .  ,  w  T  =  90°).  The 
bilinear  approximation  illustrated  (curve  i)  is  a  digitized  4th  order 
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F  I  GURE  EZ- 1  -  MOM-  RECURSIVE  FILTERS 


r 


Butterworth  filter  with  a  resulting  transfer  function 


Hr(Z) 


1  4-  4  Z'1  +  6  Z~2  +  4  Z~3  +  Z~4 
10.  64  +  5.  17  Z"2  +  .  188  Z'4 


IV -22 


The  above  equation  was  obtained  utilizing  the  transformation  technique  dis¬ 
cussed  in  Section  II.  Curves  2  and  3  are  obtained  by  applying  Eq .  IV -20  to 
Eq .  IV  -22  yielding  a  set  of  nonrecursivc  coefficients  si  vert  as 

TABLE  1 


co  = 

.  094 

C1  = 

.  376 

l2  = 

.  518 

C3 

.  193 

C4 

.  159 

o 

tl 

.  101 

C6 

.  068 

c_  = 

i 

.  045 

C8 

.  030 

C9  = 

.  020 

The  nonrecursive  filter  given  as  curve  2  utilizes  C0  through  C^.  The  filter 
of  curve  3  utilizes  hQ  through  h„  . 

Curve  4  is  obtained  utilizing  a  conventional  Fourier  series  expansion 
of  the  ideal  rectangular  filter.  Thus  the  coefficients  can  be  shown  to  be 


i 

I  where 


n 


1 

0  “  T 


h  =  h 
n  -n 


sin  n  tt  / 2 
n  tt 


H4(Z) 


u  =  0 


-u 


IV  -23 


I V  -  2  4 


i 


a 
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As  can  be  seen  from  Figure  IV  -  1  .  curves  2  and  3  represent  signifi¬ 


cantly  better  approximations  to  the  bilinear  recursive  form  than  the  conven¬ 
tional  Fourier  approximation.  The  ripple  characteristics  in  both  the  pass 
and  stopbands  have  been  reduced.  Another  comparison  of  these  various 
filter  forms  is  shown  in  Figure  IV  2.  Here  each  of  the  four  filters  accepted 
a  step  at  its  input  and  the  filter  outputs  computed.  Table  II  summarizes  the 
results . 


T  ^  PI.  E  II 


Filter  Outputs 

(for  Step  Input) 

m 

Bilinear 
(Curve  1) 

10  teem  Approx.  8  term  Approx. 

(Curve  2}  (Curve  3) 

F  ourier 
Approximation 
(Curve  4) 

y0 

.  094 

.  094 

.  094 

.  064 

yl 

.  47 

.  47 

.  47 

.  064 

y2 

.  988 

.  9«8 

.  988 

-  .042 

y3 

1.181 

1.181 

1.181 

-  .  042 

y  4 

1. 022 

1 . 022 

1 . 022 

.  276 

y5 

.  921 

.  921 

.  921 

.  776 

y6 

.  989 

.  989 

.  989 

1 . 094 

y7 

1 . 03  5 

1. 03  4 

1 . 034 

1. 094 

y8 

1. 005 

1. 004 

.  988 

y9 

.  98  4 

.  98  4 

.  988 

V,0 

.  998 

1 . 053 

^  oo 

To  1 .000 

. 
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From  this  figure,  the  improvement  in  curves  2  and  3  as  an  approximation  to 
1  (over  that  of  the  Fourier  approximation,  curve  4)  is  evident. 

The  approximation  technique  used  to  obtain  curves  2  and  3  can  be 
termed  a  "Digital  Impulse  Invariant"  technique  m  that  it  yields  an  approxima¬ 
tion  to  the  impulse  response  of  the  recursive  (bilinear  in  this  case)  filter 
which  is  exact  up  to  the  truncation  point  --  that  is,  up  to  the  last  coefficient 
retained.  This  procedure  is  the  digital  counterpart  to  the  conversirn  of  an 
analog  filter  to  a  digital  filter  through  the  use  of  standard  Z-tran  sfnrm  s , 

This  approach  yielded  a  point  for  point  match  with  the  analog  filter's  impulse 
response,  Still  another  interpretation  of  this  technique  is  that  the  bilinear 
transformation  'rounded'  the  corner  of  the  ideal  filter  of  Figure  1V-1  and 
therefore  this  impulse  invariant  technique  provides  a  simple  procedure  for 
smooth*  S3  or  modifying  the  Fourier  coefficients.  Lastly,  this  procedure  is 
simple  to  implement. 

C.  ZERO  REMOVAL  AND  RELOCATION 

As  noted  previously,  ther*1  are  often  two  sources  of  error  inherent  in 
the  design  of  nomocursive  filters.  The  second  of  these  errors  was  due  to 
the  truncation  procedure  used  to  yield  a  finite  memory  filter.  A  possible 
approach  to  nonrecursive  filter  design  is  to  view  the  problem  in  such  a  way 
that  the  second  source  of  error  appears  to  be  absent.  This  is  accomplished 
as  follows:  From  Eq.  1 V  —  1 ,  one  can  view  the  N+  1  unknown  coefficients  as 
amplitude  samples  of  so, no  time  waveform.  Thus,  this  time  waveform  is 
finite  in  extent  (all  other  samples  are  actually  equal  to  zero)  by  definition  and 
not  due  to  truncation.  Viewed  in  this  manner,  the  frequency  response 
characteristics  achievable  are  those  whose  transforms  are  finite  pulses. 
Ideally,  it  would  be  desirable  to  have  the  complete  mapping  of  all  finite  pulses 
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and  the  variety  of  frequency  responses  they  yield.  Then  one  would  choose  the 
response  "closest"  to  the  one  desired  and  then  sample  the  tinite  pulse  at  an 
appropriate  rate.  Although  the  above  mapping  was  not  accomplished ,  the 
following  will  illustrate  the  procedure  and  will  serve  to  tie  together  the 
various  approaches  to  nonreevr sive  design. 

Initially  a  characterization  of  finite  pulses  was  desired.  As  noted  by 
Campbe)',  et.al/^  a  pU]se  (in  time)  of  finite  extent  is  completely  character¬ 
ized  by  tile  fact  that  its  transform  (complex  Fourier)  has  an  infinite  sot  of 
zeros.  This  is  most  easily  seen  by  noting  a  rectangular  pulse  and  its  sin  x/x 
transform  which  consists  of  an  infinite  number  of  zeros  uniformly  distributed 
ovt-r  the  real  frequency  axis.  Subject  to  certain  constraints^  ^  these  zeros 
can  be  removed  and/or  relocated  to  alter  the  frequency  characteristic  while 
still  retaining  the  finite  pulse  structure  of  the  time  waveform. 

Perhaps  the  simplest  illustration  of  the  process  of  zero  removal  is  to 
consider  a  rectangular  pulse  existing  for  |tj  «c  1/2,  If  the  first  zero  pair  is 
removed,  it  can  be  shown  that  the  resulting  pulse  is  of  the  form 


(1  +  cos  2 -n  t)  p(t)  IV -25 

where  p(t)  is  the  rectangular  pulse  existing  over  |t|  <  !/2. 

Thus,  the  familiar  cos^x  weighting  function  can  be  viewed  in  the  fre¬ 
quency  domain  as  a  multiplication  of  the  sin  x/x  form  by  a  factor  { 1  /  1  -f  '). 
Similarly,  the  form  of  the  pulse,  as  a  result  of  removing  N  zeros  from  a 
rectangular  pulse  whose  extent  is  from  -ir  <  t  <  ti  can  be  shown  to  be 


S«>  -  -zV 


N 

V 

V  A  cos  u  r  (t  f  tt  )  ;  a  =  — —  IV -26 

—  n  m-  1  '  rt  r 

n  -  1  -> 


where 


N 

V 


IV  -27 


»  "i  >  > 

IT  r ,  r  7T  r  -  r 
k  u  k  n 


ami  the  zeros  ol  the  rectangular  pulse  arc  a  I  1  r  ,  ,  1  r  , .  •  •  •  ;  u  r  /r.  . 

1  c  Mi  n  1 

1  i  om  lY-2ii,  it  can  be  shown  that  G (o)  tin  Kiiuricr  transform  of  g(l)  is 


ink  sin  tt  (-o  -  k  ) 
l  i 


-i  n  sin  tt  {to  ^  k. ) 


r  K  1 

sin  n  (.a  -  k.) 

Ci(u>)  r  2tt  -  y  A.  e  os  it  k - 7 - r-J-  IV -28 

TT  ^  — i  1  1  TT  (.J  -  K.  ) 

1-1 


Thus,  the  shape  ol  the  pulse  in  the  frequency  domain  is  determined  by  a 
weighted  summation  of  delayed  sin  x/x  waveforms  . 

For  a  pulse  of  finite  extent,  tt  can  also  be  shown  that  its  Fourier 
transform  can  be  written  at 


C.(w)  T 


G(n)  - 


sir  tt  (-i  -  n) 


n  (->  -  n) 


IV -29 


Thus,  there  is  a  link  w  tsveen  the  amplitude  sample;  of  the  frequency  rcspor.st 
and  the  A.  coefficients  which  are  determined  by  the  zeros. 

t 

In  a  sirnilar  manner,  a  finite  cos  n  pulse  ran  be  obtained  through  the 
process  of  zero  relocation  in  the  following  manner.  Assume  a  rectangular 
pulse  c\  cr  the  interval  1 1  j  c  jt  ,  Its  transform  can  be  shown  to  be  of  the 


G  (w)  -  2  it 


IV  -30 


with  zeros  occurring  at  t  1,  t  2,  •  •  •  . 


Now,  if  the  k1^  zero-pair  is  moved  to  the  location  of  the  pair,  then 
GM  becomes 


G  j  (oj) 


sin  it 


tr  co( 


) 


The  corresponding  time  function  becomes 
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g  (t)  =  1  -  (I'1.)  cos  k{t  +  tt  )  1 1 1  <  rr  IV-32 

1  4. 


where  p  =  r/k.  Thus,  one  now  has  independent  control  over  the  amplitude  and 
frequency  of  this  time  function. 

Since  the  family  of  raised  cosine  pulses  are  used  frequently  as  weight¬ 
ing  functions,  the  transforms  of  these  pulses  were  obtained  to  note  their 
relationship  to  the  process  of  zero  removal. 

If  one  assumes  a  family  of  raised  cosine  pulses  of  the  form 

g(t)  -  COS  -r—  |t|  <  -rr 

1  "  **  IV -3  3 

=  0  otherwise 

Then  G(f),  the  Fourier  transform  of  g(t)  can  be  shown  to  be 


Off)  = 


ZT 


n'.  cos  v  f  T 


7— TJ/2' 

7 T 

k  -  0 


n  oda 


9  7  ? 

[ 2 k *r  1)"  -  4T 


IV -3  4 


or 
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G(f)  = 


Tn: 


n/2  r 

7 r 

i 


2  2  f 

4k  (1  -  T 


f  sin  tt  f  T  1 

1  '  [  tt  fT  J  ’ 


n  even 


IV -3  5 


Note,  the  form  of  I  v' -  3  b  is  such  that  a  cosine  pulse  raised  to  an  even  power  of 
n  is  equivalent  to  removing  the  first  n/2  zeros  of  sin(Tr  f  T)/tt  fT.  These 
zeros  occur  at  f  =  k/  T  for  k  =  1,  2,  •  •  •  ,  n/2  and  are  those  removed  by  the 
factor  (1  -  T"  f  t'/k<S  . 

From  F(_j.  IV -3  5,  a  link  can  be  shown  to  exist  between  the  process  of 

zero  remova1,  modified  or  weighted  Fourier  series,  and  the  technique  Dro- 

(n) 

posed  by  liltiv  kniau.  'I  his  link  is  illustrated  as  follows:  Consider  the  syn¬ 
thesis  of  an  arbitrary  pulse  given  by 

N  -  1 

f(t)  =  V  A  cosn  <  "*)  |t|  <  Jr  IV-36 

— i  n  i  —  z 

n  =  0 


This  synthesis  can  be  viewed  in  two  sections.  First  over  even  n  and  then 
ovei  odd  n.  Those  raised  cosine  pulses  over  even  n  are  equivalent  to  a 
weighted  addition  of  Zero  removed  sm  x/x  pulses,  A  cosine  pulse  raised  to 
an  odd  power  of  n  can  be  viewed  as  the  product  of  an  even  power  (n  -  1) 
multiplied  by  cos  (it  t/T).  Thus  its  transform  can  be  written  in  the  form  of 
Fq.  I V  -  3  5  with  a  shift  to  trie  left  and  right  (in  frequency)  of  (l/Z  T),  Equation 
IV-36  i therefore  a  weighted  summation  of  zero  removed  sin  (it  f  T)/ir  fT 
pulse s . 

The-  relationship  between  L.q.  IV  -  3  0  and  a  Fourier  series  expansion 
is:  obtained  by  noting  the  term  of  1  c  lie  bye  he  ft  polynomials.  These  poly¬ 


nomials  are  defined  as 


IV -37 


and  if  x  =  cos  6 


V^(x)  =  cos(n  cos 


x) 


V  (x)  =  cos  n  9  IV  -38 

n 

then  a  trigonometric  polynomial  or  Fourier  series 

N  -  1 

f{6)  =  Yj  Bn  cos  n  0  IV -3  9 

n  =  0 

is  equivalent  to  a  polynomial  in  x  given  by 

N  -  1 

f^cos  *  xj  =  ^  IV -40 

n  =  0 

The  relationship  between  the  B  and  C  coefficients  relate  the  expansion  of 

n  n 

IV-39.  That  is,  the  are  identical  to  the  for  x  -  cos(w  t)/T.  The  poly¬ 
nomials  given  by  Eq.  IV-37  relate  to  and  coefficients.  The  procedure 
would  be  to  first  expand  a  function  in  a  Fourier  series.  Then,  through  the 
use  of  the  Tchebycheff  polynomials,  calculate  the  Cr  coefficients.  Using 
equation  (or  expansion)  IV-3f>  expand  the  function  yielding  the  A  coefficients. 
Then  relate  the  and  B^  coefficients  yielding  a  weighted  Fourier  expansion. 

Blackman's  technique  involves  the  synthesis  of  a  desired  character¬ 
istic  in  the  frequency  domain  through  a  summation  of  raised  cosine  pulses  of 
the  form 

N  -  1 

I  H(cos  u»  T)  |  ^  =  V  h  cos11  wT  IV -41 

n  =  0 


where  u  is  the  digital  frequency  function.  Thus  Blackman's  technique  is  a 
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synthesis  procedure  with  raised  cosine  pulses  in  the  frequency  domain  and, 
therefore,  zero-removed  sin  x/x  pulses  in  the  time  domain.  His  procedure 
attempts  to  link  the  h^  coefficients  with  those  given  by 

N  -  1 

H(Z)  =  £  an  Z"n  IV -42 

n  =  0 

Although  the  process  of  zero  removal  and  relocation  was  not  extend¬ 
ed  to  yield  the  complete  characterization  of  finite  pulses,  it  did  illustrate 
the  tie-in  among  the  various  synthesis  techniques.  These  techniques  all 
center  about  a  finite  Fourier  series  approximation  in  some  form,  whether 
it  results  from  truncation  of  the  infinite  Fourier  series  expansion  of  the 
desired  function,  the  truncation  of  the  infinite  Fourier  expansion  of  the 
characteristic  obtained  through  the  bilinear  transform  (as  in  IV-B)  or  the 
discrete  version  discussed  in  IV -A. 

Two  basic  procedures  which  can  be  used  in  the  design  of  nonrecur¬ 
sive  digital  filters  have  been  discussed.  The  first  procedure  was  based  on 
the  assumption  that  the  desired  frequency  characteristics  were  available 
only  at  discrete  points.  The  second  approach  proceeded  on  the  assumption 
that  a  satisfactory  recursive  design  was  available.  These  recursive  coef¬ 
ficients  were  then  converted  to  a  nonrecursive  design  with  illustrative  ex¬ 
amples  indicating  the  simplicity  of  the  technique  as  well  as  improvement 
that  this  technique  offers  over  the  more  conventional  Fourier  series  approach. 
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V.  ERROR  CONSIDERATIONS  I 

I 

I 

(9)  ( i  0)  ! 

As  noted  by  Kaiser,  '  Gold'  v  and  others,  there  are  essentially  ! 

three  sources  of  errors  associated  with  digital  or  discrete  time  processing.  j 

j 

The  first  error  results  from  the  sampling  and  quantizati  .n  of  the  input  ! 

6ignal.  The  second  error  is  due  to  the  representation  of  the  filtering  coef-  | 

ficients  with  a  finite  number  of  bits.  The  third  error  is  one  of  computation-  | 

al  quantization.  That  is,  the  digital  filter  is  an  arithme'ie  unit  which  per- 

; 

forms  the  operations  of  weighted  multiplication  and  additions.  There  arc  j 

therefore  errors  caused  by  roundoff  of  these  results  which  are  often  j 

*  } 

utilized  in  further  computations.  < 

| 

In  this  section,  various  aspects  of  the  latter  two  types  of  errors  \ 

i 

t 

will  be  discussed.  The  relationship  between  coefficient  accuracy,  stabil-  ■ 

ity  and  critical  filter  parameters  will  be  detailed  as  well  as  the  effect  of 

alternate  synthesis  or  implementation  procedures.  The  relationships  i 

between  computational  quantization,  normalized  sampling  rates  and  various 
implementations  will  be  discussed  along  with  an  error  reducing  procedure. 

A.  STABILITY  AND  COEFFICIENT  ACCURACY 

It  can  be  shown  that  the  evaluation  of  the  denominator  of  the  digital 
filter  transfer  function  H(Z)  at  the  point  Z  *  =  1  gives  an  indication  of  the 
maximum  allowed  variation  of  the  digital  filter  coefficients  for  stability  to 
be  maintained.^  ^  For  a  digital  filter  of  the  form 


N 


H(Z)  = 


V  a  Z 

i~i  n 
n=  0 

T7 - 


-n 


V 

u 


b 

n 


n  =  0 


V-l 


this  index  is  of  the  form 


N 

i  =  y  b.  v-2 

u  j 
j  =  i 


or  as  it  is  more  conventional  to  divide  through  by  the  b^  coefficient,  this 
index  is  given  as 


F 


0 


N 

V 

Lj 

j=  1 


V-3 


The  tabular  approach  to  digital  lowpass  filter  synthesis  utilizing  the 
bilinear  transform  as  discussed  in  Section  II,  provides  for  a  simple  evalua¬ 
tion  and  interpretation  of  the  above  stability  index.  Consider  the  talle  used 
to  generate  the  coefficients  as  shown  on  ihe  next  page. 


TABLE  I 


Coefficients 

Pqk” 

u 

P .  K  —  ”  ^ 

1  OJ 

U 

P?K^-2 

L  to 

U 
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PN 

K° 

U 
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1 

z*1(a1-b1) 

(  ) 
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) 

-2 

z 
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) 

• 

• 

• 

d. 

1  >  J 

' 

■ 

z-N(aN’V 

(  ) 

.... 

( 

) 

where 


In  addition. 


rN  -  Ni 
r  "  rT' 

the  constraining  relationships  among  the 


d.. 

iJ 


are  given  by 


N 

T  d.  -  0  1  <  j  <  N 

‘-J  t.J  -  — 

1 


M 

V  d.  =  2*  V  -4 

CJ  1,0 
i  =  0 

To  obtain  the  denominator  (b^)  coefficients  frcm  Table  I,  one  substi 
tutes  for  Pj  the  appropriately  subscripted  denominator  coefficients  of  the 
analog  filter  (B  )  and  sums  all  the  products  of  the  table  entries  and  column 
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headings  for  row  j.  If,  however,  the  stability  index  is  formed  by  summing 


dowm  each  column  then,  utilizing  V-4,  F^  becomes 


F 


0 


<2Kn  )N  Bo 

u 

HR - 


y  b.  k 

i  =  1 


N  -i 
CJ 

u 


where 


K_ 


tanfir 


V-5 


V  -  6 


It  can  also  be  shown  that  stability  becomes  marginal  when  Fq  is  equal 

to  zero.  Thus,  the  coefficient  accuracy  problem  is  most  acute  when  the 

order  of  the  filter  is  high  and  the  factor  K_  is  small.  With  K_  defined  as 

u  u 

in  V-6,  small  values  of  K_  occur  when  the  ratio  of  sampling  rate  to  cutoff 

u 

frequency  is  high.  When  this  is  so,  V-5  can  be  approximated  by 

F„  «<2K_)N  V-7 

0  u 

u 

and  the  digital  filter  coefficients  approach  in  magnitude  the  binomial  coef¬ 
ficients  and  alternate  in  sign.  This  result  is  obtained  from  Table  I  in  that 
the  coefficients  of  the  last  column  can  be  shown  to  be  the  above  mentioned 
binomial  coefficients. 

This  equation  clearly  exhibits  the  effect  of  the  order  of  the  filter  and 
the  ratio  of  bandwidth  to  sampling  rate  on  the  stability  problem.  Table  I 
also  gives  a  great  deal  of  insight  into  the  effect  of  truncating  the  filter  coef¬ 
ficients.  Consider  row  i  of  the  table,  given  as 
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( 


R.  =  CN  B  K-  +  d  B  K^'  1  +  -  •  •  d  D  KIl"a  +  1  + 
i  i  0  to  i  ,  1  w  i ,  a  a  -  1  u> 

u  u  \i 


t  cN  bn 

1  N 


Now  the  coefficients  cf  the  analog  filter  paired  as  (By,  B^.),  (Bj,  B^  _  j) 

<B.,  BN  -  .)  are  most  often  of  the  same  order  of  magnitude,  Furthermore,  it 

can  be  shown  that  the  d.  ,  coefficients  can  also  be  paired  in  the  same  manner. 

1  >  J 

Thus,  if  the  summation  is  accomplished  in  the  foregoing  pair-wise  manner, 
then  one  is  adding  terms  of  the  form 


M  K  —  KN.'2ktl 

u  u> 

u  u 


When  the  ratio  of  sampling  rate  to  bandwidth  is  high  and  the  oriiir  ot  the 

filter,  (N)  is  large,  then  K^-  ^  can  be  considerably  less  than  1  --  especially 

u 

for  small  k  (towards  the  outer  edges  of  the  table).  Therefore,  in  effect,  the 
truncation  problem  effects  the  highest  order  analog  filter  coefficient  first. 
That  is,  Bn  is  effectively  set  equal  to  zero  in  both  the  numerator  and  de¬ 
nominator.  Thus,  if  the  analog  filter  was  originally  of  the  form 


G(S)  - 


ansN  +  an-isN‘I  +  A( 


BmS,n  +  B, 
N  I 


V  -  1  0 


It  then  becomes 


Ci(H)  - 


N  -  1 

AN  -  1  S  +  '  '  '  Ao 


BN  -  1  s 


'^Although  the  numerator  and  denominator  are  usually  oi  ditierent 
order,  they  can  always  be  written  as  shown,  by  adding  the  required  zeri 


Therefore ,  one  is  actually  synthesizing  a  lower  order  filter  as  one  truncates 
the  coefficients.  Of  course,  the  more  coarse  the  quantization,  the  more 
coefficients  (from  the  high  end)  become  deleted.  This  reduction  of  the  order 
of  the  filter  explains  the  results  obtained,  in  that  as  the  quantization  process 
continues,  the  filter  slopes  do  decrease,  which  signifies  a  lower  order 
filter.  The  onset  of  instability  is  due  to  the  fact  that  the  remaining  coef¬ 
ficients  are  not  changed  as  the  order  is  reduced  so  that  the  zeros  (or  poles) 
of  the  polynomial  shift  in  a  complicated  manner.  It  should  be  noted  that  for 

values  of  K_  close  to  .  5,  Fn  gets  smaller  with  increasing  N  at  a  very  slow 
u 

rate.  This  value  of  K__  corresponds  to  a  sampling  to  cutoff  frequency  ratio 

wu 

of  approximately  7:1. 

A  similar  result  can  be  obtained  for  the  bandpass  case.  This  is 

most  easily  seen  if  one  applies  the  LP-BP  transformation  with  a  4:1  ratio  of 

sampling  to  carrier  frequency.  As  noted  in  Section  II,  this  transformation 

- 1  -2 

is  obtained  by  replacing  Z  by  -Z  and  K_  by  Kg  .  Thus,  the  same  table 

u 

can  be  used  if  alternate  rows  (even  numbered)  are  multiplied  by  -1  and 

K_  — *•  K-5  .  It  can  also  be  shown  that  the  sum  of  the  entries  in  any  partic- 
w  13 

u 

ular  column  is  equal  to  zero  except  for  column  1  where  the  entries  sum  to 
N 

2  .  Since,  in  the  low-pass  case  stability  is  determined  by  evaluating  the 
denominator  of  the  transfer  function  at  Z  ^  =  1,  then  the  low-pass  case 
yields  Bn  2N  K  —  .  For  the  bandpass  case,  the  point  Z  =  e  replaces 

U  (x) 

u 

the  point  Z  *  =  1  in  the  stability  computations.  However,  for  the  case 

1 

under  discussion,  w  T  =  tt/2  and  Z  =  "£  j.  Thus,  Z  =  -1  for  n  odd  and 
- 

2  '  =  1  for  n  even.  This  alternation  in  sign  effectively  cancels  out  the 

alternation  in  sign  due  to  the  transformation  and  the  stability  index  for  band¬ 
pass  filters  is  given  by  By  2^  .  Thus  the  index  becomes 


S3 


In 

V 

Li 


„  ,.N  -  i 

BiKii 


l  -  0 


where  the  B.  coefficients  are  the  denominator  coe flieient s  of  the  low-pass 
analog  filter. 

The  foregoing  results  indicate  that  synthesis  of  a  high  order  filter  in 
conjunction  with  a  high  ratio  of  sampling  to  cutoff  (or  bandwidth)  frequency 

should  be  avoided.  These  conditions  on  the  factor  K  would,  for  example, 

oy 

vi 

occur  whenever  it  is  desired  to  filter  a  narrow  band  segment  from  a  broad 
bandwidth  signal.  The  narrow  band  constraint  dictates  a  small  whereas 
the  broad  bandwidth  dictates  (through  the  use  of  the  sampling  theorem)  a 
large  f  .  The  foregoing  results  indicate  that  rapid  degradation  occurs  when 
the  above  ratio  exceeds  7:1  (K ...  <  .5).  Under  these  conditions,  a  cascade 

VX>  - - 

u 

approach  to  digital  filter  synthesis  yields  a  significant  improvement.  The 
following  will  illustrate  this  improvement. 


B.  CASCADE  vs.  DIRECT  SYNTHESIS 

With  the  foregoing  constraints  to  be  illustrated,  an  elliptic  filter 
was  designed  with  the  specifications  assumed  iu  he 

Pa  s  sband  Spec .  .01  db  ripple  0<co<i 

Stopband  Spec.  40db  attenuation  1.38  <  co  <  oo 

Utilizing  the  techniques  of  Section  III,  tlie  above  specifications  led  to  a  tixMi 
order  analog  filter  of  the  form 


8  4 


A  Z 

S  +  A,S  +  A. 
t-  u 

6 

V  B.s1 


V-  1 3 


where  the  coefficients  art  given  as 


TABLE  II 


Aq  =  6.  7  9609  B^  =  5.7  17  37 


Az  -  5.  40108  B  ^  =  6.  29689 


Bn  ^  1.26743  B  =  5.  06663 

vJ  4 


13  j  -  3.44569  B5  =  2.  62193 


B6=  1 


Solving  the  numerator  and  denominator  polynomials  of  V-13  the  pole -zero 
configuration  for  this  filter  can  be  shown  to  be 


TABLE  HI 


zQ  =  t  j  1 . 84538 
z  =  t  j  1  41268 


.  10102  1  j  1.  12369 
. 39699  t  j  .9921  1 


(zeros) 


(poles) 


p3  =  -  . 81294  t  j  . 45946 


Initially,  this  filter  was  digitized  under  the  conditions  that  the  signal  ’S 
sampled  at  6  times  the  cutoff  frequency  ol  the  filter.  Thus,  K  =  tan(u>  3T/2) 
tan(ir/6)  =  .  5773  5.  Utilizing  Table  li  of  Section  II,  substitution  of  the  above 
parameters  yields  tnc  following  digital  filter  coefficients. 


\m  IllHIWIli  r  -i:.: 


TABLE  IV 


Numerator  C oe  ft . 

a  -  ]  .  1851  589 

?  t  =  2.  0438!  1 5 

a 2  =  2.  8421442 

a.  -  3. 9669833 
o 

a.  =  2. 8421442 
4 

a5  =  2. 04381 15 
afe  =  1  .  1851589 


Denominator  Coeft. 

b0  =  6.  3  1772532 

b,  =  -12. 9964865 
1 

bz  -  18. 4185255 

b3  =  -1  4.  8467182 

b  .  -  8. 34145098 

4 

b  =  -  2. 65466173 
b6  -  . 424433844 


This  filter  and  the  filter  obtained  by  rounding  the  above  coefficients  to  five 
significant  figures,  three  to  the  right  and  ivvo  to  the  left,  are  plotted  in 
Figure  V-l.  As  dictated  by  the  6:1  ratio,  the  cutoff  frequency  of  this  filter 
is  at  2r/6rad  or  60°.  As  can  be  seen,  there  is  essentially  no  difference 
in  the  filter  responses. 

This  filter  was  then  synthesized  using  the  cascade  approach.  The 
procedure  used  was  to  perform  the  factorization  in  the  analog  domain  in  the 
foil  owing  manner:  '1  he  .liter  G(S)  was  expressed  in  cascade  form  as 


G(S)  r  G,(S)  G  2  ib)  G  ,  (b) 


-  1  4 


where  both  G^(S)  and  G2<8>)  arc  of  the  form 


C*.j(S)(G,(.S)) 


r  A 

li. S  i  '■'i 


V  •  1  5 


The  filter  to..  65)  is  of  the  form 

j  ' 


T&UKlCATtO  COEFF, 


1 


V-  1  6 


G3<S) 


S*  + 


BjS  + 


B, 


Using  the  results  of  Table  III,  the  foregoing  filter  characteristics  become 


^  /Ck  _  <S  +  j  1 . 84538)(S  -  j  1. 84538) 

=  Iff +  TTTnD2'+  j  1.  1  i* 3 69T(S" +~1  OTCZ- -  ' j  1 !  12T69) 


V-  17 


r  ,cv  _  (S  +  j  1  ■  41 2 68) (S  -  j  1.41268) _ 

C2(b)  "  (5  +  .  39699  +  j  .  9921  0(S  +  •  396M  -  j  .  99211) 


V-  18 


C3(S)  =  (5  +  .K1294  +  j  .  4B94S)(5  -T !  HT2~94  -  j  .  45TO  V_1 

If  Eqs .  V-17  through  V-19  are  expanded  and  the  bilinear  transform  applied 

to  these  filters  utilizing  the  technique  of  Section  II-A  with  K_  =  .57735  = 

u 

tan(3  0°)  then  the  three  digital  filters  become 

Filter  1  (from  Gj(S)) 


a0 

=  2. 

135 

ii 

o 

1. 

541 

al 

=  . 

27  0 

bl  = 

-1 . 

151 

a2 

=  2. 

135 

b2 

1. 

308 

Fil 

ter  2  (from 

<S)) 

a0 

=  1. 

665 

b0  = 

1 . 

839 

al 

= 

67  0 

bl  = 

-1 . 

239 

a2 

=  1. 

,  665 

b2 

922 

88 


Filter  3  (from  G^(S)) 

aQ  =  .  333  bQ  =  2. 230 

=  .  667  =  -1.419 

a2  =  . 333  b2  =  . 352 

The  product  of  these  responses  was  obtained  and  compared  with  that  of 
Figure  II- 1,  showing  no  significant  difference. 

The  above  procedure  was  repeated  for  a  higher  ratio  of  sampling  to 
cutoff  frequency  to  illustrate  the  previously  mentioned  degradation.  Using 
K_  =  .  17633  =  tan(10°)  (i.e.  ,  a  20°  cutoff)  the  sixth  order  digital  filter  coef 

OJ 

u 

ficients  become 

TABLE  V 


aQ  =  .036517931 

a  =  -.050516111 

a2  =  -.033249521 

a3  =  . 107569042 

a  =  033249521 

a.  =  -.050516111 
o 

a6  =  .02  651793  1 


bQ  =  1.66053236 

bj  *  ■  8.  1507299 

b2  =  17. 0485004 

b3  =  -19.  3912146 

b  =  12.  626522 

4 

bc  =  -  4. 4568362 

D 

b6  =  .665664076 


This  filter  and  the  one  obtained  by  rounding  the  coefficients  as  before,  (3 
decimal  digits)  are  plotted  in  Figures  V-2  and  V-3.  Note  the  severe  change 
in  filter  performance.  That  is,  the  curve  of  V-3  is  unrecognizable  as  a  low 
pass  digital  filter.  This  illustrates  the  problems  involved  in  a  direct 
synthesis  of  a  high  order  filter  where  the  ratio  of  sampling  rate  to  cutoff 
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frequency  is  large  and  coefficient  rounding  is  desired.  Once  again  this 
filter  was  synthesized  as  three  cascaded  filters  with  the  respective  digital 
coefficients  obtained  as  before  as  shown  below. 


Filter  1 

(from  G  j 

(S)) 

=  1.106 

bo 

1. 

175 

=  -1.788 

bl  = 

-1 . 

921 

1! 

O 

O' 

b2 

1. 

004 

Filter  2 

(from 

(S)) 

=  1.062 

b0  = 

1. 

175 

=  -1.876 

bl  = 

-1. 

929 

=  1.062 

b2 

895 

Filter  3 

(from  G^ 

(S)> 

=  .031 

bo  = 

1. 

314 

=  .  062 

bl* 

-1 . 

946 

=  .  031 

b2  = 

74U 

?  was  obtained  and 

plotted 

■  and  that 

of  the  rounded  a 

advantage 

of  the  ca 

iscade  a 

An  alternate  approach  to  comparing  the  cascade  vs.  direct  approach 

( 121 

was  investigated  by  Knowles  and  Oicayto,  '  who  modeled  the  coefficient 
truncation  problem  as  a  parallel  error  filter  and  used  a  statistical  approach 
to  evaluate  performance. 
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The  foregoing  results  were  expected  in  view  of  the  results  on  coef¬ 
ficient  errors  and  generalized  design  procedures  which  showed  that  the 

N  -  j 

coefficients  are  primarily  determined  by  factors  R.K _  where  R,  is  the 

1  co  j 

J  u  J 

analog  filter  coefficient,  N  is  the  order  of  the  filter  and  is  the  normal- 

u 

ized  cutoff  frequency.  Thus,  as  I<_  is  reduced  the  effect  of  R.  is  lost  for 

1  to  ) 

u  J 

large  N  as  the  coefficients  are  rounded. 

C.  COMPUTATIONAL  QUANTIZATION 

The  nature  of  computational  quantization  is  dependent  upon  the  partic¬ 
ular  implementation  used.  A  digital  filter  recursion  equation  expressed  in 
direct  form  is  given  as 


N 

y  =  ^  a .  x 

n  -j  j  n  -  j 


N 

V  b.  y 


J  n  -  j 


V  -20 


j=0  J=1 

When  computational  quantization  occurs  at  each  operation  of  V-16,  the 
recursion  equation  becomes 


N 

-  — 

N 

F  I 

=  Qm 

v  Q'  . 

a .  x 

-  V  Q" 

b  y 

V-21 

11 

—  n,  j 

J  n  -  j 

— '  n,  j 

J  n-j 

j  =  o 

- 

j=  > 

L  J 

where  ana  Q"[*3  represent  the  operation  of  quantizing  the 

function  on  which  it  operates. 

If  one  represents  the  operation  of  quantization  by  an  equivalent  addi¬ 
tive  noise  source  n'  as 


Q'  [a]  =  n1  +  a 


then  V-21  becomes 


V  -22 


N  N 

y  =  Q  +  }  a.  x  V  by  , 

7  n  n  u  jn-j  1—  j  7  n  -  1 


j  =  0 


j=  1 


where  Qr  is  the  total  quantization  noise  in  the  system  (i.e.,  =  n'  +  n"H  n1"). 


Thus,  the  error  becomes 


y  -  y  =  c  =  Q 
n  n  n  n 


N 

^  b  e 

i  n  -  i 

J--  1  ' 


V-23 


If  it  is  assumed  that  the  quantization  noise  has  a  white  spectral  density,  the 
mean  square  output  noise  is  given  by 


M  .  S .  O .  N 


<rrr>$i>Vz>K,,r'1>¥ 


V  -24 


where 


HjfZ)  - 


i  -  y  b  z 

~ i  i 


V  -  25 


i  =  1 


2  . 


the  b^  are  the  denominator  coefficients  ol  the  original  filter  and  Oq  is  the 
mean  square  input  noise  which  is  dependent  upon  the  degree  of  quantization. 

Thus,  computational  quantization  has  been  expressed  by  an  additive 
noise  source,  processed  by  H^(7.)  which  consists  only  ol  the  denominator  cr 
poles  of  the  original  filter. 

Consider  the  input-output  equations  ol  a  digital  filter  synthesized 
using  the  canonical  form.  In  this  case  the  tr  .  tfer  tunaion  corr esponding 
to  V-20  given  as 
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y  a  z'n 

u  n 


H,7v  -  Y<Z)  -  n  =  0 

H(Z)  -  TT2 y  “ rr- 


-  N<z> 

■  ’em 


y  b  z'n 

Z/  n 


is  rewritten  as 


V  -26 


X(Z)  _  Y (  Z)  _ 

m  -  nITT  F(Z> 


V  -  27 


where  F(Z)  is  some  arbitrary  function.  Inverse  transforming  the  two  equa¬ 
tions  specified  in  V -27  yields 


f  =  x  -  y  b.  f 
n  n  4_i  j  n  -j 

j=  1 

V-28 


y  =  y  a .  f 
n  <-<  J  n  -  j 


where  f  is  an  intermediate  variable  introduced.  This  is  the  so-called 
n 

canonical  form. 

If,  as  before,  one  assumes  that  quantization  occurs  at  each  stage  and 
that  the  operation  of  quantization  is  represented  by  an  additive  noise  source, 
then  it  can  he  shown  that  the  error  is  given  by 


N  N 

£  -  y  a  Q  .  -  ^  b  £ 

n  — t  1  a  -  j  J  n  -  j 

j  =  0  j  1 


V  -29 


This  pair  ol  equations  is  related  to  the  state -space  concept  utilized 


in  ontima)  control. 
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Thus,  under  a  white  noise  assumption,  the  mean  square  output  noise  is 
given  by 

2. 

M.S.O.N.  =  {-2VJ  )  H2(Z)  H^Z'1)  -1^  V  -30 

where 


Thus,  for  this  form,  the  noise  is  processed  by  II0(Z)  which  is  identical 

£- 

to  the  original  filter  H(Z).  Hj(Z)  consisted  of  only  the  poles  of  the  original 
filter  whereas  H,(Z)  consists  of  the  poles  and  zeros,  that  is  the  entire  original 
filter . 

Utilizing  the  relationships  concerning  the  generalized  design  proced¬ 
ures  for  lowpass  filters  (Table  I)  it  was  desired  to  determine  the  relationship 
between  the  output  noise  and  the  ratio  of  sampling  rate  to  cutoff  frequency 
for  lowpass  filters  synthesized  using  both  of  the  previous  approaches. 

The  filters  investigated  were  second  order  lowpass,  all  pole  filters 
of  the  form 

G(S)  =  - - - - -  V  -32 

AS  +  BS  +  C 

Evaluation  of  the  integrals  (M .  S . O .  N . /a^  )  and  utilization  of  the  generalized 
synthesis  techniques  yielded 
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Direct  Synthesis 


C.  Ki.  +  A 

M.S.O.N.  _  wu 

- 2 - - - T 

o~  16ABCK- 

Q  oj 


V  -33 


Canonical  Synthesis 


where 


M.S.O.N. 


~~T 

aQ 


K  _  (B  K—  +  A) 

U  OJ 

_ u _ u _ 

CkI  +  BK_  +  A 

u  oj 

u  u 


V  -3  4 


V  -3  5 


The  above  results  were  evaluated  for  a  Butterworth,  Chtbycheff  ()/2db 
ripple)  and  a  Maximally  Flat  Time  Delay  Filter.  The  coefficients  for  these 
filters  are 


Butterworth 


A  =  1  B  - JT 


P] 


V-36 


Chebycheff 

A  =  1  B=  1.4256245  C=  1.5162026 


p  7128122  +  j  1 .0040425  V-37 

Maximally  Flat  Time  Delay 

A  =  1  B  =  3  C  =  3 

Pj  =  -  1 . 5  +  j . 8660254 


V  -  38 


These  results  are  plotted  in  Figures  V-4,  5  and  6. 

It  is  to  be  noted  that  at  high  ratios  of  sampling  rate  to  cutoff  frequency, 

the  canonical  form  is  vastly  superior.  There  are,  however,  regions  where 

the  direct  synthesis  should  be  used, 

The  results  obtained  indicate  the  importance  of  determining  the  value 

of  the  parameters  K  „  .  This  parameter  not  only  determines  the  synthesis 

wu 

technique  used,  but  also  determines  the  required  coefficient  accuracy, 

D.  AN  AUXILIARY  STORAGE  T ECH1QUE 

In  the  previous  section,  expressions  were  derived  showing  the  effect 
of  computational  quantization  as  an  error  producing  equivalent  noise  source. 
This  discussion  will  develop  the  concept  of  an  auxiliary  storage  to  be  used 
to  reduce  these  errors. 

If  *he  round-off  portion  discussed  previously  is  stored  in  an  auxiliary 
storage  and  the  main  computation  proceeds  as  prescribed  by  the  applicable 
model  and  some  quantization  levei,  the  accuracy  of  the  output  value,  y  ,  may¬ 
be  corrected  by  appropriately  weighting  the  contents  of  the  storage  and  add¬ 
ing  it  to  the  output. 

In  order  to  make  a  "perfect"  correction,  all  prior  round-off -seg  - 
merits  have  to  be  stored  and  multiplied  by  appropriate  coefficients  with 
complete  pre  ;ision.  Since  this  is  obviously  impractical,  in  the  general 
rase,  the  technique  developed  will  permit  a  conservative  estimate  to  be 
made  of  the  error  due  to  truncation  in  previously  rounded-off  samples  as 
well  as  the  quantization  present  in  the  round-off  arithmetic. 

From  £q.  V-2.',  the  error,  E(Z),  and  the  round-off  or  quantization 


9d 


Q!  Z)  are  related  by 


FIGU  RE  3Z-5:  MAXIMALLY  FLAT  HUE'  DELAY 


loo 


E(Z)  = 


Q( z) 


N 

1  +  l  bj  z'j 

j=  1 


V  -39 


As  noted  previously,  this  recursive  equation  can  be  converted  to  a  non- 
r  ecursive  form  yielding 


Thus 


where 


V 

£  =  ;  u-,  q  . 

n  <-j  k  n  -  k 

k  =  0 


oo 

y  =  y  +  T*  C.  q  . 

7  n  7  n  L- 1  k  nn  -  k 

k  =  0 


C0=  1 


N 

C.  =  -  V  b.  C.  . 
k  Li  j  k  -  j 

j=  1 


V  -40 


V  «  4 1 


V-42 


A  block  diagram  of  the  system  utilizing  Eq.  V-40  and  truncating 
the  nonrecursive  form  after  k  terms,  is  illustrated  in  Figure  V  -  7  . 

There  are  three  ways  to  use  this  technique  to  advantage,  In  the  first 
instance,  where  only  the  filter  output  is  necessary,  a  relatively  gross 
running  arithmetic  would  be  necessary,  after  which  the  ptn  output  would  be 
refined  using  some  number  of  current  and  prior  round-off  values  taken  from 
the  piggy-bank.  This  simplifies  the  running  arithmetic  at  the  expense  of 
increased  storage. 


lOf 


L 
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F\C?  URE  31-7:  A  R£CROS\V£  DIGITAL  F-U.TE.R  WITH  AM  AUXILIARY  STORAGE. 


The  second  use  of  the  storage  is  to  periodically  test  the  error,  (at  a 
slower  than  input  sample  rate)  improving  the  value  y^  with  the  weighted  over¬ 
flow  whenever  the  error  value  overflows.  Overflow  is  defined  as  that  part 
of  the  error  equal  to  or  greater  than  the  gross  quantization  of  the  running 
arithmetic.  The  residue  of  the  overflow  would  be  reinserted  as  the  new 
round-off  for  the  n  value  and  the  whole  process  repeated. 

The  third  possible  application  is  dependent  on  the  nature  of  the  over¬ 
flow.  If  the  bandwidth  of  the  overflow  samples  is  small  relative  to  the  input 
samples,  then  it  is  conceivable  that  a  refinement  of  y^  can  be  made  periodi¬ 
cally  but  at  a  rate  slower  than  the  input  sampling  rate. 

An  Illustrative  Example 

An  example  is  next  presented  for  an  arbitrary  filter  where  the  sum 
never  overflows.  However,  a  correction  is  made  at  the  81*1  sample  time 
using  a  truncated  (approximate)  and  non-truncated  (exact)  technique. 

Example 

Assume  the  filter  is  given  by 

y  anx  -b.y  -b-y  , 
n  On  1  n  -  1  2  n  -  2 

where : 

ao  1 

bj  =  1/2 

b2  =  1/4 

x  =  1  n  >  0 

n  — 

x  =  0  n  <  0 

n 

Table  VI  lists  exact  valueB  of  the  output,  y  ,  approximate  values  of  the  out¬ 
put,  y  ,  assuming  a  quantization  of  l/8  and  the  resulting  round-off. 
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TABLE  VI 


EXACT  AND  APPROXIMATE  OUTPUTS  VALUES 


n 

Qnx 

0  n 

Vn-l 

b2yn-2 

yn 

yn 

Vn-l 

Vn-2 

% 

0 

i 

0 

0 

1 

1 

0 

0 

0 

1 

i 

1/2 

0 

1/2 

»/2 

1/2 

0 

0 

2 

i 

1/4 

1/4 

1/2 

1/2 

1/4 

1/4 

0 

3 

i 

1/4 

1/8 

5/8 

5/8 

1/4 

1/8 

0 

4 

i 

5/16 

1/8 

9/16 

l/2 

5/16 

1/8 

1/16 

5 

i 

9/32 

5/32 

9/16 

1/2 

1/4 

5/32 

3/32 

6 

i 

9/32 

9/64 

37/64 

5/8 

1/4 

1/8 

0 

7 

i 

37/128 

9/64 

73/128 

1/2 

5/16 

1/8 

1/16 

8 

i 

73/256 

37/256 

73/128 

1/2 

1/4 

5/32 

3/32 

9 

i 

73/512 

73/512 

293/512 

5/8 

1/4 

1/8 

0 

Assume  a  quantization  of  1  / 8 
and  record  (i.e.,  store)  the 
round-off  q 

_ jn _ 

By  using  the  algorithm  for  obtaining  C  ,  the  following  coefficients 
are  developed. 

CQ  =  1,  CL  =  1/2,  C2  =  0,  C3  =  -1/8,  C4  -  -1/16,  C6  =  0,  C6  =  l/64, 

C  =  1/128,  •• •  . 

Using  these  coefficients,  the  error  for  the  8th  output  will  be  computed. 
E„  =  !  x  3/32  -  1/2  x  1/16  +  0x0+  l/8  x  3/32  -  l/lfe  x  1  / 1 6  +0x0 

O 

+  1/64  x  0  •  • 

E0  =  24/256  -  8/256  +  3/256  -  1/256  =  18/256  =  9/126 

o 


y8  =  1/2  =  64/128 

since 

yg =  ys +  Es 

y8=  73/128 

The  above  numerical  example  illusirate6  that  the  exact  result  is  ob¬ 
tained  if  all  of  the  round-off  values  are  used. 

If  only  the  current  and  next  to  last  round-off  values  are  stored,  an 
approximate  correction  to  y^  is  obtained  as  y^ 

where  y  =  y  +  E 
7  n  ’  n  n 

Eq  =  1  x  3/2  -  1/2  x  1/16  =  2/32  =  8/128 
o 

y8  =  72/128. 

The  error,  after  refinement,  using  the  approximate  correction  is  one 
part  in  128  while  the  error  between  yg  and  yg  is  nine  parts  in  128. 

If  the  round-off  error  is  pessimistically  assumed  to  be  always  equal 
to  the  least  significant  bit  (LSB)  then  the  max  error  max  is  defined  as 

E  (max)  =  lim  S(LSB) 
n 

k  — ♦  ou 

where  „ 

lim  S  =  ^  C 

^  g 

k  — ♦  ao  g  =  0 

For  the  example  above  lim  S  =  4/7  and  LSB  -  1/8. 

k  — ♦  oo 

.  '  .  max  4/7  x  l/8  =  1  / 1  4 .  If  the  storage  values  are  not  utilized 
at  all,  the  max  error  for  the  conditions  assumed  would  be  less  than  l/l4. 
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If  the  stored  round-off  error  itself  were  quantized  (which  is  generally 
the  case)  then  the  error  (conservative)  due  to  this  approximation  would  be 

lim  S  •  (LSB).  Where  (LSB)  is  the  least  significant  bit  of  stored  round-off 
k  — >  oo 

value.  Thus,  this  auxiliary  storage  technique  offers  the  possibility  of  reduc¬ 
ing  errors  . 

Two  sources  of  errors  in  digital  filtering  have  been  discussed.  The 
errors  resulting  from  a  truncation  of  the  filter  coefficients  were  related  to 
such  critical  filter  parameters  as  the  order  of  the  filter  and  the  ratio  of 
sampling  rate  to  critical  frequencies  using  the  tabular  approach.  This  ap¬ 
proach  allowed  for  an  interpretation  of  the  filter  degradation,  the  onset  of 
instability  and  the  bit  truncation  of  the  filter  coefficients.  The  conditions 
under  which  the  cascade  synthesis  approach  is  preferred  is  also  discussed, 
The  errors  due  to  computational  quantization  and  the  effect  of  both  the  direct 
and  canonical  implementation  forms  on  these  errors  are  then  discussed. 
Curves  are  obtained  showing  the  regions  in  w'hich  one  implementation  is 
superior  to  the  other.  Lastly,  an  error  reduction  technique  utilizing  an 
auxiliary  storage  approach  is  presented.  Inis  tec  unique  offers  the  possibili¬ 
ty  for  reducing  computational  quantization  errors. 
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VI.  RECOMMENDATIONS 


As  is  often  the  case  with  analytical  studies,  the  various  techniques 
investigated  suggest  areas  where  further  efforts  would  be  desirable.  In 
addition,  some  of  the  results  themselves  may  show  sufficient  promise  to 
warrant  experimental  justification  or  breadboarding.  Summarized  below 
are  the  areas  which,  as  a  result  of  this  investigation,  fall  into  the  above 
categories . 

1.  The  results  of  Section  II  (B  and  C)  and  Section  III,  concerning  the 
various  techniques  and  constraints  which  are  associated  with  the 
design  of  bandpass  filters,  point  out  the  interplay  between  the  samp¬ 
ling  rate,  carrier  frequency,  and  filter  bandwidth.  The  foregoing, 
in  addition  to  the  results  of  Appendix  A,  clearly  shows  the  need  to 
develop  sampling  and  processing  techniques  which  will  enable  one  to 
operate  digital  filters  at  sampling  rates  at  or  near  the  information 
bandwidth. 


The  results  of  Section  V  (A  and  C)  concerning  stability  and  coef¬ 
ficient  representation  illustrate  the  difficulties  encountered  when 
filters  are  to  be  designe',  under  the  constraints  of  small  K_  factors 

and  high  order  filters.  Jmce  truncation  invites  stability  problems 
and  the  foregoing  constraints  require  an  inordinate  amount  of  pro¬ 
cessing,  techniques  should  be  investigated  which  overcome  tl.ese 
difficulties  and  will  thus  allow  for  the  ocessing  or  extraction  of 
narrow  hand  information  from  broadband  signals. 


3.  The  error  considerations  discussed  in  Section  III  indicate  that  the 

development  of  techniques  which  provide  for  approximate  filter  shapes 
through  coefficient  modification  should  be  investigated.  These 
modifications  should  be  directed  towards  an  overall  reduction  in  the 
time  required  to  perform  the  arithmetic  processing. 


4.  The  results  of  Section  II  (B  and  C)  and  those  of  Section  III  (C.  and  K) 
suggest  that  a  more  detailed  study  be  made  of  the  relative  advantages 
of  the  two  specific  bandpass  filter  design  techniques  discussed.  The 
shifting  technique  offers  a  possibility  for  developing  the  bandpass 
coefficients  from  those  of  the  lowpass  filter  directly  on  line.  It  also 
appears  more  desirable  for  tracking  filter  applications.  The  LP-BP 
transformation  o;  rs  other  advantages  in  terms  of  control  over  band¬ 
width  and  carrie-  requency. 


5.  The  results  of  Section  IV  indicate  that  a  further  investigation  be  direct¬ 
ed  toward  the  relationship  between  the  process  of  zero-removal  and 
relocation  and  the  class  of  allowable  frequency  response  characteristics 
achievable  for  nonrecursive  filters. 


6.  An  extension  of  the  results  of  Section  V -C  concerning  computational 
quantization  is  desirable.  This  extension  to  higher  order  filters  as 
well  as  arbitrary  bandpass  filters  would  allow  the  designer  to  proper¬ 
ly  choose  the  best  implementation  for  the  particular  critical  para¬ 
meters  of  the  filter. 

7.  The  error  reduction  technique  discussed  in  Section  V-D  has  shown 
sufficient  merit  to  warrant  an  experimental  br eadboar ding .  This 
auxiliary  storage  technique  area  requires  some  additional  study  with 
regard  to  the  optimum  number  of  coefficients  to  be  used  in  the  storage 
"bank"  as  well  as  a  determination  of  which  of  the*  three  inodes  it  should 
operate. 

8.  An  experimental  investigation  should  be  initiated  concerning  the  digital 
oscillator  implementation  technique  discussed  in  Appendix  B.  This 
technique  has  applications  to  such  areas  as  frequency  translators, 
bandpass  processing  and  tracking  filters. 

9.  Besides  the  discussion  of  the  shifting  technique  for  bandpass  digital 
filters  and  the  nonrecursive  design  from  tabulated  data,  the  majority 
of  the  synthesis  techniques  discussed  have  utilized  a  transformation 
procedure  which  converts  an  analog  filter  to  an  equivalent  digital 
filter.  These  approaches  might  be  generally  termed  digital  equivalence 
procedures.  An  alternate  approach  is  to  synthesize  digital  filters  with¬ 
out  reference  to  an  equivalent  analog  filter.  Thi3  approach,  which 
might  be  termed  direct  digital  synthesis,  warrants  detailed  investiga¬ 
tion  in  that  it  will  allow  »or  more  flexibility  in  the  designs  that  can  be 
achieved.  Two  approaches  to  this  problem  of  direct  design  which  have 
been  investigated  to  a  limited  extent  under  very  specific  conditions 
have  utilized  polynomial  approximation  procedures  and  pole -zero 
shifting  procedures  in  the  Z  plane.  Roth  of  these  procedures  as  wel.‘ 

as  several  other  approaches  should  be  expanded  for  use  in  the  direct 
synthesis  of  digital  filters. 
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APPENDIX  A 


A  BANDPASS  SAMPLING  TECHNIQUE 

The  sampling  and  reconstruction  of  a  lowpass  signal  with  a  (positive) 
bandwidth  of  B/2(H?.)  and  no  spectral  energy  beyond  this  frequency  can  be 
accomplished  (theoretically)  by  sampling  at  a  rate  of  B  samples  per  second. 

B  also  represents  twice  the  highest  frequency  component  present  in  the  signal. 
If  chis  procedure  were  carried  over  directly  to  bandpass  signals,  unreasonably 
high  sampling  rates  would  result.  Moreover,  since  it  is  well  known  that  the 
information  content  in  a  signal  is  dependent  upon  its  bandwidth  and  not  its 
center  frequency,  one  would  expect  that  it  is  possible  to  recover  the  bandpass 
waveforem  with  a  sampling  procedure  at  rates  in  the  order  of  the  bandwidth 
of  the  bandpass  signal.  Direct  application  by  the  sampling  theory  for  band¬ 
pass  signals  show’s  this  to  be  the  case. 

The  periodic  nature  of  digital  filters  suggests  the  following  sampling 
procedure.  Assume  a  bandpass  signal  as  shown  below. 


i 

*  C. 


»  i 


-A/Vv 


MAG, 

A 


i 


1  K- 


The  process  of  sampling  produces  spectral  repeats  or  aliased  spectra 

at  multiples  of  the  sampling  frequency.  The  repeats  are  shown  as  the  dotted 

spectra.  It  is  desired  to  choose  a  sampling  rate,  fs,  such  that  there  are  no 

spectral  overlaps.  With  this  objective  the  following  constraints  can  be  placed 

til 

upon  fs,  B,  fj  and  f^.  Let  k  be  the  index  representing  the  k  repeat  or 
shifted  spectrum.  Then  for  no  overlap  at  bandpass  one  requires 


-  f,  +  k  f  >  f,  ; 
2  s  2 

Simplifying  equation  A-l  yields 


-  fj +  £s  <ft 


£ 

s 


f 


2 


2B 
k-  1 


or 


f 

s 


IT 


f  7  U  2 

s  .  ^  2  L 

T  iTl  T  T7i 


A-l 


A-2 


A-  3 


From  these  equations  one  obtains 


f  >  2B 

s  A-4 

^2 

x<  ^  <  4 

Thus,  the  absolute  minimum  allowable  sampling  rate  is  2B.  The 
allowable  normalized  sampling  rates  relative  to  f2/B»  the  ratio  of  the  upper 
frequency  to  the  bandwidth,  are  shown  in  figure  A-l.  This  diagram  shows 
that  care  must  be  exercized  in  choosing  a  sampling  rate.  It  is  of  interest  to 
note  that  the  sampling  rate  of  2B  is  allowed  only  when  i'2  and  B  are  integrally 
related.  It  is  also  of  interest  to  note  that  if  two  sampling  rates  fg  j/B  and 
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{  - /B are  chosen 
s  2 ' 

shown  to  be 


(f^,  >  £^),  then  the  allowable  sampling 


2B 


2 

. 


region  or  area  can  be 


A-  5 


which  is  independent  of  k.  Furthermore,  the  extend  of  validity  on  the  f^ /B 
axis  is 


—  - 

‘f  -  f 

f 

'2 

k 

s  2  si 

+ 

"eT 

=  1 

.  B 

2B 

which  is  an  increasing  function  of  k. 

The  foregoing  considered  the  required  •  ampling  rate  for  no  bandpass 
overlap.  Additional  constraints  should  be  placed  on  the  sampling  rate  in  order 
to  filter  this  signal  digitally  at  a  reasonably  low  frequency.  One  such  con¬ 
straint  is  to  require  that  there  be  no  spectral  interchange.  That  is,  if  the 
repeat  is  the  last  negative  repeat  from  the  negative  spectral  lobe  (moving  to 

f 

the  right)  the  (r+l)“  repeat  is  the  last  positive  repeat  from  the  positive  spec¬ 
tral  h  be  (moving  to  the  left);  then  for  no  spectral  interchange  one  requires 
that 

(r  +  1)  f  -  f,  >  f,  -  r  f 
s  2  2  s 

A-7 

f  ?  f, 

W  >'ZTTT  "B" 


As  an  example  of  the  case  where  the  sampling  rate  is  sufficient  for  no  spectral 
overlap  but  spectral  interchange  does  occur,  consider  the  following.  A  band¬ 
pass  signal  has  f  =  250  Hz,  fj  =  150  Hz,  B  =  100  as  shown  below 


1  13 


7  ( HZ) 

f  ia  chosen  as  270  Hz.  As  shown  in  the  figure,  no  spectral  overlap  occurs 
but  spectral  interchange  occurs  since 


0  1 
*•100-4  UlOO-J 


^  =2,70 


=  2.  50 


r  =  0 


A-f 


Thus  from  A-7,  2.  7  must  be  greater  than  2x2.  5  =  5  which  is  not  the  case. 


APPENDIX  B 


DIGITAL  OSCILLATORS 

Digit'll  oscillators  plsy  an  integral  part  in  such  areas  as  frequency 
translators,  heterodyning  techniques  and  tracking  and  frequency  hopping 
filters  as  well  as  other  applications  .  Thus  several  configurations  for  the 
digital  generation  of  sine  and  cosine  oscillators  will  be  discussed. 

A  starting  point  in  the  development  of  these  generators  is  their  Z 
transforms  which  can  be  shown  to  be 


ros  ui  t 
o 


Z  I  cos  u  t 
o 


sin  u:  t 

o 


Z  f  sin  u  t 
o 


1  -  (cos  a)  Z 


-1 


1  -  2{cos  a)  Z'1  +  Z‘* 


(sin  a )  Z 


-1 


TT 

1  -  2<cos  a)  Z  +  Z 


~2 


B-l 


B-2 


where  a  is  2rr  divided  by  the  ratio  of  sampling  frequency  to  center  frequency. 
That  is  the  sinusoidal  values  are  obtained  at  increments  of  a  (radians). 

In  equations  B-l  and  B-2  the  Z  transforms  are  in  the  form  of  a  retio. 
Therefore,  they  can  be  interpreted  as  transfer  functions  of  a  recursive 
digital  filter,  inverse  trans forming  these  equations  leads,  respectively,  to 

cos  u  t 
o 

yk  =  *k  '  <COS  a)  XK-1  +  (2  CO*  a)yk-l  -  yk-2 


Sill  ii)  t 
O 


yk  -  (sin  a)  xk_,  t  (2  Cos  a)  y^,  -  Vk  2 


1  5 


B-4 


Thus,  these  generators  have  been  modeled  as  recursive  digital  filters  whose 
impulse  responses  provide  the  samples  of  sine  and  cosine  at  any  desired  in¬ 
crement.  That  is,  if  20  samples  per  cycle  were  desired,  then  a  =  (2tt/20)  =  18° 
and  the  outputs  would  correspond  io  the  values  of  these  sinusoids  at  18°  incre¬ 
ments.  It  should  also  be  noted  that  these  digital  filters  have  poles  on  the  unit 
circle  and  thus  if  allowed  to  run  indefinitely  an  "  infinite"  amount  of  noise 
would  result  due  to  computational  quantization.  This  problem  can  be  circum¬ 
vented  by  periodically  re-starting  the  oscillators.  The  implementation  of  these 
generators  is  shown  in  figure  B-l.  The  dotted  region  of  the  figure  is  shown 
merely  for  illustrative  purposes  and  represents  the  first  two  initial  conditions 
to  be  loaded  into  the  delay  elements,  That  is,  the  cosine  generator  requires 
K  and  K  cos  a  and  the  sine  generator  0  and  K  sin  a .  Thus  B-3  and  B-4  can  both 
be  written  as 

cos  oj  t  (and  sin  w  t) 

_ o _ ' _ o  _ 

yk  =  (2  cos  a)  yk  >  *  ^  B-5 

where  the  cosine  terms  are  obtained  by  setting  (y  y  ., )  to  be  (K(  K  cos  a)  and 
the  sine  terms  by  setting  the  initial  outputs  to  be  (0,  K  sin  a). 

The  configuration  shown  in  figure  B-l  is  instrumented  in  direct  form 
and,  therefore,  has  some  of  the  error  problems  discussed  in  section  VI  with 
regard  to  computational  quantization.  An  alternate  configuration  is  based  on 
the  identities 

cos  (A-t-  B)  -  cos  A  cos  B  -  sin  A  sin  B 
sin  (At  B)  =  sin  A  cos  B  +  cos  A  sm  B 
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B-  6 


OUTPUT 


(b'  SINE  GE NER^TOR 
P  iOURE  6- t 


These  identities  lead  to  the  recursion  equations 


y  -  (cos  q)  y  .  +  (sin  a)  x  , 

'  n  n-1  n-1 

B  -  7 

x  =  (cos  a)  x  .  -  (sin  a)  y  , 
n  n-1  ■'n-1 


where  the  y*  s  represent  the  sine  output  and  the  x*  s  the  cosine  output.  This 
configuration  is  shown  in  figure  B-2.  For  sampling  rates  that  are  in  the  order 
of  10-30:1  and  resetting  the  system  every  cycle,  these  two  configurations  pro¬ 
duced  very  similar  results.  Without  resetting,  the  configuration  of  figure  B-2 
was  superior.  A  considerable  simplification  occurs  in  both  of  these  configura¬ 
tions  when  these  generators  yield  4  samples  per  cycle.  In  this  case, 
a  =  n/2  =  90°,  cos  u  =  0  and  sin  a  =  1  and  the  outputs  can  be  shown  to  yield 
the  sequence  1,  0,  -1,  0,  1,  0,  -1,  ....  This  of  course  merely'  corresponds  to 
a  sampled  square  wave  and  there  are  no  computational  quantization  problems. 
For  this  case,  the  oscillators  have  been  reduced  to  a  rather  trivial  configura¬ 
tion.  This  result  is  expected  in  that  it  is  common  practice  to  build  hetero- 
dyner9whieh  square  wave  modulate  a  signal  in  order  to  shift,  the  signal  by  an 
amount  (in  the  frequency  domain)  equal  to  the  fundamental  frequency  of  the 
square  wave.  Energy'  centered  at  the  higher  harmonics  is  then  filtered  out. 

A  third  configuration  yielding  improved  performance  with  less  hardware 
is  shown  in  figure  B-3.  This  configuration  is  based  on  expressing  13-1  and  B-2 
in  canonical  form.  The  resulting  equations  are 


f  -  x  +  (2  cos  a)  i  .  -  f 
n  r  n  - 1  r  -  2 


••  =  (sin  a)  :  . 

n  n  -  i 


f’  =  x'  H  (2  cos  a)  f  .  -  f 
n  n-i  n  -  2 

=  1 ,,  -  ( ->s  o  )  f]% 

1  IS 


B-t 


From  the  above  it  can  be  seen  that  since  x'  =  x  ,  f '  =  f  that  the  sine 

n  n  n  n 

and  cosine  generators  sh?.re  the  same  auxiliary  variable  but  are  not  interlocked 
as  in  figure  B-2.  It  is  also  to  be  noted  that  there  is  a  reduction  in  the  required 
number  of  multipliers  using  sine  a  and  cos  a.  Furthermore,  the  errors  in  the 
sine  a  multiplication  do  not  affect  the  alternate  or  cosine  output. 

Using  the  techniques  discussed  in  section  V,  it  was  shown  that  the 
error  performance  of  this  configuration  is  superior  to  that  of  the  previous 
configurations.  Simulation  of  the  recursion  equations  with  10  bit  and  14  bit 
quantization  bears  this  out  as  shown  in  tables  I  and  II. 

TABLE  I 

10  Bit  Quantization 

Angle  (degrees)  Generated  (sine)  Ideal  (rounded) 

(Figure  B-2)  (Figure  B-3) 


0  0  0  0 


10 

.  174 

.  174 

.  174 

20 

.  343 

.  342 

.  342 

30 

.  5C1 

.  501 

.  500 

40 

.  644 

.  644 

.  643 

50 

.  768 

.  767 

.  766 

60 

.  869 

.  867 

.  866 

70 

.  948 

.  942 

.  940 

80 

.  989 

.  987 

.  985 

90 

1.  006 

1.  003 

1.  000 

TABLE  II 


14  Bit  Quantization 


Angle  (degrees)  Generated  (cosine)  Ideal  (rounded) 


0 

1.  000 

1.  0000 

10 

.  9848 

.  9848 

20 

.  9401 

.  9397 

30 

.8671 

.  8660 

40 

.  76^9 

.  7660 

50 

.  6458 

.  6428 

60 

.  5043 

.  5000 

70 

.  3478 

.  3420 

80 

.  1807 

.  1737 

90 

.  0083 

.  0000 

The  results  shown  in  this  table  are  for  output  generated  at  10°  incre¬ 
ments.  For  smaller  increment,  the  improvement  shown  in  the  performance 
of  the  configuration  of  figure  B-3  over  that  of  the  previous  figures  increases. 
This  is  once  again  due  to  the  fact  that  the  canonical  form  is  best  when  the 
sampling  rate  is  high.  The  outputs  shown  in  these  tables  result  from  quan¬ 
tizing  the  appropriate  equations  to  the  specified  number  of  bits. 

Digital  oscillators  not  only  are  important  devices  as  frequency  syn¬ 
thesizers  but  also  play  an  integral  part  in  such  important  communication 
devices  as  translators,  tracking  filters  and  bandpass  processing  among  many 
others.  The  configuration  of  figure  B-3  can  therefore  be  combined  with  the 
results  concerning  digital  bandpass  filters  to  yield  a  completely  digital  device. 
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